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COSMIC  RAY  EFFECTS  ON  MICROELECTRONICS 


Part  III:  Propagation  of  Cosmic  Rays 
in  the  Atmosphere 


1.0.  Introduction: 

This  report  considers  the  propagation  of  primary  cosmic  rays  through 
the  atmosphere  and  assesses  their  effectiveness  at  producing  single  event 
upsets  ( SEU )  in  microelectronic  components.  The  problem  of  soft  upsets 
centers  on  the  ability  of  a  single  charged  particle,  upon  passage  through 
a  sensitive  volume  on  a  silicon  chip,  to  cause  the  logical  state  of  a 
memory  circuit  to  change.  Such  upsets  obviously  have  immense  leverage 
In,  for  example,  control  circuit  applications.  Shielding  against 
cosmic-ray-induced  upsets  is  difficult  because  they  are  highly 
penetrating  and  have  measurable  fluxes  at  extremely  high  energies.  The 
upset  potential  of  cosmic  rays  Is  dependent  on  many  factors  including  the 
solar  cycle,  geomagnetic  field  strength,  atmospheric  depth,  and  device 
orientation.  In  this  report  upset  calculations  in  the  atmosphere  are 
described  and  computer  code  to  perform  these  calculations  Is  presented. 

In  Section  2  we  consider  the  effect  of  passage  through  air  on  cosmic 
ray  fluxes.  The  most  prominent  interactions  with  air  are  nuclear 
fragmentation  reactions  and  ionization  loss.  Close  collisions  between 
air  nuclei  and  cosmic  rays  cause  breakup  of  the  cosmic  ray  into  lighter 
nuclides.  Generally  a  few  nucleons  or  alpha  particles  are  emitted.  The 
fragments  retain  roughly  the  same  energy/nucleon  as  the  incident 
nucleus.  Secondaries,  that  is,  accelerated  fragments  of  silicon  nuclei 
are  not  considered  in  this  report  although  they  are  a  potential  source  of 
upsets  (Peterson,  1981).  Leptons  and  mesons  created  In  the  collisions 
are  also  not  considered. 

Distant  collisions  between  air  molecules  and  cosmic  rays  do  not 
change  the  composition.  They  generally  cause  ionization  of  the  air. 
These  collisons  gradually  decelerate  the  cosmic  rays  and  modify  their 
spectrum.  The  rate  of  energy  loss,  or  stopping  power,  is  proportional  to 
the  square  of  the  cosmic  ray  charge.  Thus  the  heavy  ions  which  are  most 
effective  at  causing  upsets  suffer  greatest  attenuation  In  the 
atmosphere.  Nuclear  decay  need  not  be  considered  because  propagation 
through  the  atmosphere  takes  a  fraction  of  a  millisecond  while  most 
decays  occur  on  time  scales  of  more  than  milliseconds. 

In  Section  3  results  of  cosmic  ray  propagation  through  air  are 
Integrated  to  give  total  flux  levels  at  any  altitude.  This  integration 
involves  calculating  flux  levels  as  a  function  of  zenith  and  azimuthal 
angles.  These  depend  on  the  geomagnetic  field  and  the  amount  of  air  to 
the  top  of  the  atmosphere  in  any  direction.  The  latter  is  a  sec  e  effect 
for  the  altitudes  considered  in  this  report.  The  geomagnetic  field 
effect  results  in  an  enhancement  of  fluxes  to  the  west  relative  to  the 
vertical.  Total  cosmic  ray  fluxes  as  a  function  of  altitude  and  cutoff 
are  displayed  showing  the  roughly  exponential  decrease  as  a  function  of 
both  variables. 

Section  4  considers  the  conversion  of  cosmic  ray  spectra  into  LET 
(linear  energy  transfer)  spectra  and  upset  computations.  LET  spectra 
show  the  flux  as  a  function  of  stopping  power  in  silicon.  Instead  of 
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classifying  according  to  energy,  they  are  classified  according  to  the 
rate  at  which  they  give  up  energy  In  silicon.  As  stated  previously  LET 
increases  as  the  square  of  the  cosmic  ray  charge.  LET  Is  not  a  monotonic 
function  of  energy.  In  fact  the  lowest  LET  occurs  at  high  energies  (a 
few  GeV/nucleon)  while  the  highest  LET  occurs  at  a  few  MeV/nucleon.  The 
stopping  powers  vary  by  roughly  a  factor  of  50  in  this  range.  LET 
spectra  for  geomagnetic  cutoffs  ranging  from  0  to  8  GV,  and  altitudes 
between  55,000  and  150,000  ft.  are  presented. 

Given  the  device  parameters,  LET  spectra  may  be  converted  into  upset 
probabilities.  In  any  environment  the  probability  of  upset  depends  on  a 
critical  charge  and  sensitive  volume  (device  dependent  parameters)  In 
addition  to  the  cosmic  ray  LET  soectrum.  The  critical  charge  Is  the 
number  of  free  electron-hole  pairs  necessary  In  the  neighborhood  of  a 
junction  to  cause  an  upset.  These  pairs  are  created  at  a  rate  of  one  per 
3.6  eV  of  deposited  energy.  The  sensitive  volume  is  characterized  by  a 
chord  length  distribution  giving  the  relative  frequency  of  each 
pathlength.  The  energy  deposited  is  proportional  to  this  pathlength. 

Section  5  consists  of  descriptions  of  the  programs  developed  to 
estimate  soft  upsets  in  the  atmosphere. 


2.0.  Propagation  of  Cosmic  Rays  through  Air 
2.1  Master  Propagation  Equation 

The  cosmic  ray  flux  at  the  top  of  the  atmosphere  consists  roughly  of 
2500  proton/ (m2  sec  ster),  about  one  tenth  as  many  He  nuclei,  and 
roughly  1  per  cent  heavier  nuclei.  The  average  energy  of  these  particles 
is  1  or  2  GeV/nucleon.  Within  the  magnetosphere  and  atmosphere  they  are 
subject  to  several  attenuating  processes.  The  most  Important  are: 

1)  Fragmentation  in  collision  with  N  and  0  nuclei  in  the  air. 

2)  Energy  loss  by  Ionization  of  the  air. 

3)  Deflection  from  the  atmosphere  by  the  geomagnetic  field. 

The  first  two  of  these  processes  are  discussed  in  this  section,  while  the 
third  is  discussed  In  Section  3. 

The  master  equation  for  cosmic  ray  propagation  with  Ionization  loss 
and  fragmentation  (Ginzburg  and  Syrovatskl 1 ,  1964)  is: 

TT  *  “  nVo1J1+  nv.£1  °1jJj  +  h  ('dT  *7)  (2.1-1) 

Terms  In  this  equation  are  defined  as  follows: 

Ji  -  Differential  energy  spectrum  of  cosmic  rays  of  species  1 

In  part1cles/m2  sec  steradian  (MeV/nucleon).  If 
species  are  ordered  by  Increasing  mass,  and  within  that 
by  decreasing  charge,  then  all  reactions  reduce  the 
species  Index.  For  example,  ®li  <  7ge  <  ?L1. 
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oi  »  Total  reaction  cross-section  for  nuclei  of  species  I  In  air. 

a i j  •=  Partial  fragmentation  cross-sections  for  nuclei  of  species  1 
to  be  produced  in  collisions  of  species  j  with  air. 

n  ■  Number  density  of  N  and  0  nuclei  In  air. 

v  -  Cosmic  ray  velocity  which  is  related  to  energy/nucleon 

(MeV/nucleon)  by 

(WO2.  ?JJLi_1863> 

(E  ♦  031. S)2 

A-j  a  Atomic  mass  of  species  i. 

It  is  convenient  to  re-express  eq.  2.1-1  in  terms  of  pathlength, 

Instead  of  time.  The  pathlength  is  the  total  number  of  grams  of  material 
In  a  cylinder  of  1  cm2  cross  section  traversed  in  time,  t,  by  the 

particle.  Thus, 

x  (g/cm2)  »  n  v  m  t  (2.1-2) 

where  m  is  the  average  mass  of  an  air  atom.  The  pathlength  between  the 
top  of  the  atmosphere  and  any  point  directly  below  is  essentially  equal 
to  the  air  pressure  at  that  point.  However,  pathlength  has  units  of  mass 
per  unit  area  while  pressure  is  force  per  unit  area.  The  relation 
between  pathlength  and  altitude  is  shown  in  Fig.  2.1.  Data  for  this 

figure  were  taken  from  the  1962  U.  S.  Standard  Atmosphere  (NASA  et  al., 
1962).  A  fit  to  the  data  valid  between  50  and  150  kilofeet  Is  also 
shown.  Table  2.1  shows  altitude  vs.  pathlength  comparisons  of  special 
Interest  In  this  report.  The  composition  of  the  atmosphere  is 
essentially  constant  In  the  region  of  interest.  We  assume  a  mixture 

consisting  of  79  per  cent  N  and  21  per  cent  0  by  number. 

Table  2.1 

ALTITUDE  vs.  PATHLENGTH 
(corresponding  to  computations  in  this  report) 


Altitude  Pathlength 


kilofeet 

kilometers 

g/cm2 

150 

45.7 

1.33 

125 

38.1 

3.59 

100 

30.5 

11.1 

85 

25.9 

22.0 

75 

22.9 

35.0 

65 

19.8 

56.5 

60 

18.3 

72.1 

55 

16.3 

92.2 
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In  terms  of  pathlength  Eq.  2.1-1  may  be  rewritten: 


3J.|  o.|  J.j 

~W  “  “  23354 


♦  E 


j>i 


gjj_  _Jj 
23954 


(3.1-3) 


where  w-|  is  the  stopping  power  of  nuclei  of  species;  'n  air. 

This  fixes  the  cross  section  units  to  millibarns  and  the  stopping 
power  units  to  MeV/(g/cm2).  The  factor  23954  is  the  average  mass  in 
grams  of  an  air  nucleus  divided  by  the  unit  conversion  from  cm2  to 
mb. 


For  the  purposes  of  further  manipulation  we  express  Eq.  2.1-3  in 
the  condensed,  matrix  form: 

3J  tii  +  3  ( w  J)  (2.1-4) 

•Sx  "  M  J  +  IF 

The  fragmentation  effects  are  contained  in  the  modification  matrix,  M, 
and  ionization  loss  effects  in  the  derivative  term;  w  is  the  energy 
less  per  nucleon. 

Approximating  the  Ionization  loss  and  backsubstituting  for  K  we  have 
the  approximate  solution:  .  _ 

^  U  t  I  "d 

J(E,x+ax)  -  eMixj  [l  -  w(E)  J(E.x)  +  w  (E  +  AE)-||  J(E  +  aE,x)j 


In  this  form  the  flux  at  a  given  pathlength  and  energy  depends  only  on 
the  flux  at  shorter  pathlengths  and  greater  (or  equal)  energies. 

The  flux  at  the  highest  energy,  where  ionization  loss  can  be 
linearized  (Letaw  et  al.,  1983b),  is: 

J(Emax»  x)  *  e^x  J(Emax»  0)  '  (2.1-6) 

where  M  now  contains  the  term: 


-(na 


-  n  +  1) 


A1Rp 


(2.1-7) 


Rp  is  the  proton  range  in  air,  Zi  the  charge  of  species  i,  n  the 
range-energy  index  (~  1),  and  a  the  spectral  index  (~  2.7)  of  the 
differential  flux  at  high  energy.  The  relative  importance  of  this  term 
is  gauged  by  comparison  with  fragmentation  losses. 
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2.2  List  of  Isotopic  Species 

The  most  abundant  cosmic  rays  have  charges  between  1  and  28.  Heavier 
elements  exist  in  much  smaller  quantities  because  they  are  beyond  the 
stability  point  at  iron  for  stellar  nucleosynthesis.  A  list  of  230 
Isotopes  from  H  through  Fe  was  extracted  from  the  Table  of  Isotopes 
(Lederer  and  Shirley,  1978).  Most  of  these  isotopes  are  unstable  but 
cannot  decay  during  the  fraction  of  a  millisecond  spent  in  the  atmosphere. 

To  reduce  the  length  of  the  list,  a  propagation  through  air 
(neglecting  ionization  loss)  with  pathlength  of  30  g/crn^  was 
performed.  Elemental  abundances  were  taken  from  Adams,  et  al.  (1981). 
Isotopic  abundances  were  estimated  from  recent  data  (Webber,  1981;  Young 
et  al.,  1981;  Garcia-Munoz  et  al.,  1981).  Unstable  isotopes  are  created 
in  fragmentation  reactions.  Those  which  did  not  achieve  an  abundance 
greater  than  or  equal  to  1  percent  that  of  Fe  were  dropped  from  the 
list.  Those  which  developed  abundances  between  1  percent  and  5  percent 
of  iron  were  regrouped  with  other  Isotopes  of  the  same  charge.  In  this 
way  a  reduced  list  of  104  isotopes  was  developed  (Table  2.2). 


2.3  Cross  Sections 

Two  sets  of  cross  sections,  the  total  reaction  cross  sections  cr-j 
and  the  partial  fragmentation  cross  sections  <rjj,  are  needed  to  perform 
air  propagation  calculations.  These  cross  sections  are  for  reactions  of 
the  form: 


Nuclide  +  <  W  anything 


Roughly  5500  partial  cross  sections  and  104  total  cross  sections  are 
required.  These  cross  sections  are  Incompletely  known  and  are  estimated 
using  semi-empirical  formulas. 

The  partial  cross  sections  for  many  proton-nucleus  reactions  are 
given  by  the  semi-empirical  formulas  (Siiberberg  and  Tsao,  1973a, b; 
Silberberg  and  Tsao,  1977a;  Tsao,  et  al.,  1983).  These  are  scaled  up  to 
air-nucleus  interactions  using  procedures  discussed  in  Silberberg  and 
Tsao  (1977b).  The  scale  factor  is  roughly  2  with  additional  enhancement 
for  light  products. 

The  total  cross  sections  are  calculated  using  the  Bradt-Peters  (1950) 
overl ap  model : 

V  -  49.88  (At  !/3  +  Ap  1/3  -  0.4)2  (2.3-1) 

where  At  and  Ap  are  the  target  and  projectile  masses  respectively. 
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Table  2.2 


Isotope 

ProRani 

.tors 

Fraction  of 
Elemental  Abundance 

1H 

1 

5B 

9 

1.000 

1H 

2 

0.000 

2HE 

3 

0.000 

1H 

3 

0.000 

2HE 

4 

4BE 

8 

4BE 

a 

5B 

9 

5B  9 

1  .  000 

3LI 

5 

0.000 

2HE 

5 

0.000 

3LI 

6 

0.  500 

2HE 

6 

0.000 

4BE 

7 

0.  500 

3LI 

7 

3LI 

8 

3LI 

9 

0.  500 

2HE 

8 

0.000 

4BE 

9 

0.  500 

5B 

10 

5B 

8 

0.2  50 

4BE 

10 

4BE 

11 

4BE 

12 

0.000 

6C 

11 

6C 

9 

6C 

10 

0.000 

5B 

11 

5B 

12 

5B 

13 

5B 

14 

0.750 

6C 

12 

1.000 

6C 

13 

6C 

14 

6C 

15 

GC 

16 

0.000 

7N 

14 

7N 

12 

7N 

13 

0.  500 

80 

15 

80 

13 

80 

14 

0.000 

7N 

15 

7N 

16 

7N 

17 

7N 

18 

0. 500 

80 

16 

80 

17 

80 

18 

80 

19 

80  20 

1.000 

9F 

17 

0.000 

9F 

18 

0.000 

1  ONE 

19 

1  ONE 

18 

1  ONE 

17 

0.000 

9F 

19 

1.000 

10NE 

20 

0.  500 

9F 

20 

0.000 

11NA 

21 

11NA 

20 

0.000 

10NE 

21 

0.2  50 

9F 

21 

9F 

22 

0.000 

11NA 

22 

0.000 

10NE 

22 

10NE 

23 

10NE 

24 

0.2  50 

12MG 

23 

12MG 

22 

0.000 

11NA 

23 

1.000 

12MG 

24 

0. 500 

11NA 

24 

■ 

o.ooo 

13AL 

25 

13AL 

24 

0.000 

12MG 

25 

0.250 

11NA 

25 

11NA 

26 

0.000 

13AL 

26 

0.000 

12MG 

26 

12MG 

27 

1  2MG 

28 

0.2  50 

14S I 

27 

14SI 

26 

0.000 

13AL 

27 

1.000 

14SI 

28 

1.000 

13AL 

28 

0.000 

14SI 

29 

0.000 

13AL 

29 

1  3AL 

30 

0.000 

15P 

30 

ISP 

28 

1  5P 

29 

0.  000 

14SX 

30 

14SI 

31 

14SI 

32 

0.000 

15P 

31 

1.000 
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Table  2.2  (Cont’d) 


165  32  16S  30  16S  31  1.000 

15P  32  0.000 

16S  33  0.000 

15P  33  15P  34  0.000 

17CL  34  17 CL  33  0.000 

16S  34  0.000 

17CL  3  S  0. 500 

16S  35  16S  36  1 6S  37  0.000 

18AR  36  18AR  34  18AR  35  0.500 

17 CL  36  0.000 

18AR  37  0.000 

17CL  37  17CL  38  17CL  39  0.500 

19K  38  19K  37  0.000 

18AR  38  0. 500 

19K  39  0.250 

18AR  39  0.000 

2 OCA  40  2 OCA  39  0.500 

19K  40  0.250 

18AR  40  18AR  41  0.000 

2 OCA  41  0.000 

19K  41  0.500 

2 OCA  42  0.250 

19K  42  19K  43  0.000 

21SC  43  21SC  42  0.000 

2 OCA  43  0.000 

21SC  44  0.000 

2 OCA  44  2 OCA  45  20CA  46  0.250 

22TI  45  23V  45  22TI  44  0.000 

21SC  45  1.000 

22TI  46  0.  500 

21SC  46  0.000 

22TI  47  0.250 

21SC  47  2 ISC  48  0.000 

23V  48  23V  46  23V  47  0.000 

22TI  48  0.2  50 

23V  49  0.250 

22TX  49  2  2TI  50  0.000 

24CR  50  24CR  48  24CR  49  0.000 

23V  50  0.2  50 

24CR  51  0.000 

23V  51  23V  52  23V  53  0. 500 

2 5MN  52  2  5MN  51  0.000 

24CR  52  0.  500 

26FE  53  0.000 

2  5MN  53  0.2  50 

24CR  53  0.2  50 

2 6FE  54  0.000 

2  5MN  54  0.2  50 

24CR  54  0.2  50 

2 6FE  55  0.000 

25MN  55  0.500 

26FE  56  1.000 
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(2.3-2) 


When  the  projectile  Is  hydrogen  the  cross  section 

0  =1  44.9  aO-7 

(Letaw  et  al.,  1983a)  has  been  used. 

Significant  unknowns  in  the  cross  section  data  base  are  cross 
sections  for  production  of  light  products  (A  <  6).  These  cross  sections 
are  not  known  even  for  proton-nucleus  collisions.  Errors  of  20  percent 
or  more  in  the  proton  and  He  abundances  may  result  from  ignoring  the 
production  of  these  nuclides  in  fragmentation  reactions.  Rough  estimates 
of  these  cross  sections  for  production  of  light  nuclei  were  made  based  on 
data  obtained  by  Freie*  and  Waddington  (1975). 

We  note  here  that  uncertainties  in  the  total  cross  sections  are  the 
single  greatest  uncertainty  in  these  calculations.  The  error 
(**  10  percent)  in  the  total  cross  section  for  Iron  fragmentation  causes 
at  least  n  times  that  error  in  the  flux  after  n  mean  free  paths  (1  mean 
free  path  -  16  gm/cm2).  The  error  in  partial  cross  sections,  when 
averaged  over  the  isotopes  is  similar  in  magnitude  to  that  of  the  total 
cross  sections. 

We  have  neglected  the  variation  of  cross  sections  with  energy  in  this 
computation  and  used  the  high-energy  values.  This  does  not  result  in 
significant  error  because  the  mean-free -path  for  fragmentation  at  low 
energies  is  much  smaller  than  the  range. 


2.4  Ionization  Loss 


Loss  of  energy  in  the  atmosphere  due  to  ionization  of 
is  represented  by  the  stopping  power  of  the  nuclide.  The 


air  molecules 
stopping  power 


w 


dE 

Hx 


(2.4-1) 


is  the  rate  of  energy  loss  [MeV/(g/cm2)]  of  material  traversed.  For 
protons  the  stopping  power  in  air  varies  from  roughly  200  MeV/ (g/cm2) 
at  1  MeV  down  to  about  2  MeV/(g/cin2)  at  high  energies.  For  other 
particles  it  scales  as  Z2.  Table  2.3  contains  stopping  powers  of  H, 
He,  C,  0,  Ar,  and  Fe  in  air.  This  table  was  produced  by  programs 
discussed  in  Adams,  et  al .  (1983a). 


The  range  of  particles  in  air  is  also  required  in  this  work.  Range 
is  defined  by  the  integral: 


(2.4-2) 


and  is  calculated  in  g/cm2.  Table  2.4  contains  ranges  of  H,  He,  C,  0, 
Ar  and  Fe  in  air.  Ranges  for  other  nuclides  may  be  interpolated  from 
this  table. 
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Table  2.3 


Stopping  Power 


in  Air  (MeV/(gm/cm 


2 


)) 


MeV/u 

H 

He 

C 

0 

Ar 

Fe 

1.0 

2 . 205E+02 

8 . 433E+02 

5.  514E+03 

8 . 479E+03 

2. 378E+04 

3 . 472E+04 

1.2 

1.9  51E+02 

7. 502E+02 

5.  159E+03 

8. 034E+03 

2 . 343E+04 

3. 484E+04 

1.  5 

1. 673E+02 

6. 501E+02 

4. 703E+03 

7 . 437E+03 

2.277E+04 

3 . 461E+04 

2.0 

1. 363E+02 

5. 394E+02 

4. 103E+03 

6. 610E+03 

2 . 1 54E+04 

3 . 372E+04 

2.  5 

1.157E+02 

4. 652E+02 

3 . 641E+03 

5. 9  50E+03 

2 . 032E+04 

3 . 2  56E+04 

3.0 

1 . 008E+02 

4. 108E+02 

3. 275E+03 

5.  411E+03 

1 .919E+04 

3. 134E+04 

4.0 

8 . 046E+01 

3. 348E+02 

2 . 727E+03 

4. 584E+03 

.1 . 723E+04 

2 . 900E+04 

5.  0 

6.8  50E+01 

2.885E+02 

2. 361E+03 

4. 013E+03 

1 . 567  E  +  04 

2 . 69  5E+04 

6.0 

5.93  5E+01 

2. 524E+02 

2. 074E+03 

3.  559E+03 

1 . 43  6E  +  04 

2 . 513E+04 

7.0 

5.  2  52E+01 

2. 249E+02 

1 . 8  51E+03 

3. 200E+03 

1.325E+04 

2. 355E+04 

8.0 

4. 721E+01 

2 . 031E+02 

1. 673E+03 

2. 908E+03 

1 .232E+04 

2.215E+04 

10.0 

3 . 948E+01 

1 . 703E+02 

1. 406E+03 

2. 464E+03 

1 .081E+04 

1 . 982E+04 

12.0 

3 . 411E+01 

1 . 465E+02 

1 . 218E+03 

2. 142E+03 

9.647E+03 

1.796E+04 

15.0 

2 . 856E+01 

1.206E+02 

1. 020E+03 

1. 799E+03 

8 . 330E+03 

1 . 577E+04 

20.0 

2.277E+01 

9 . 301E+01 

8.  130E+02 

1.436E+03 

6.81 3E+03 

1  .  317E+04 

25.0 

1 . 908E+01 

7. 659E+01 

6. 833E+02 

1. 20BE+03 

5. 799E+03 

1 .135E+04 

30.0 

1 . 647E+01 

6. 594E+01 

5. 921E+02 

1. O49E+03 

5. 081E+03 

1 . 003E+04 

40.0 

1. 308E+01 

5.  233E+01 

4. 707E+02 

8. 3  56E+02 

4. 123E+03 

8. 249E+03 

50.0 

1.095E+01 

4. 382E+01 

3. 944E+02 

7. 005E+02 

3 . 490E+03 

7 . 077E+0  3 

60. 0 

9 . 486E+00 

3. 796E+01 

3. 418E+02 

6. 073E+02 

3 . 041E+03 

6. 218E+03 

70.0 

8 . 416E+00 

3 . 368E+01 

3. 033E+02 

5. 391E+02 

2. 708E+03 

5. 564E+03 

80.0 

7. 597E+00 

3. 040E+01 

2. 738E+02 

4.B6BE+02 

2.450E+03 

5.  053E+03 

100.  0 

6. 424E+00 

2. 571E+01 

2. 316E+02 

4.  119E+02 

2.079E+03 

4. 307E+03 

120.0 

5. 622E+00 

2.250E+01 

2.028E+02 

3. 606E+02 

1.824E+03 

3 . 7B9E+03 

150.  0 

4. 804E+00 

1 .923E+01 

1.733E+02 

3. 083E+02 

1. 562E+03 

3. 254E+03 

200.  0 

3. 967E+00 

1. 588E+01 

1 . 432E+02 

2.  547E+02 

1 . 293E+03 

2.700E+03 

2  50. 0 

3.456E+00 

1  .  383E+01 

1. 247E+02 

2. 220E+02 

1 . 128E+03 

2. 359E+03 

300.  0 

3. 113E+00 

1 . 246E+01 

1 . 124E+02 

2. OOOE+02 

1 . 017E+03 

2. 129E+03 

400.  0 

2. 684E+00 

1 . 074E+01 

9. 692E+01 

1. 725E+02 

8. 791E+02 

1.840E+03 

500.0 

2 . 430E+00 

9.726E+00 

8. 777E+01 

1.  562E+02 

7.958E+02 

1. 668E+03 

600.0 

2. 265E+00 

9 . 067E+Q0 

8. 182E+01 

1. 457E+02 

7 . 422E+Q2 

1.  557E+03 

700.0 

2. 152E+00 

8. 612E+00 

7. 773E+01 

1.  384E+02 

7 . 0  53E+02 

1 . 480E+03 

800.0 

2. 070E+00 

8 . 286E+00 

7. 479E+01 

1.  331E+02 

6.788E+02 

1 . 424E+03 

1000.0 

1. 965E+00 

7 , 864E+00 

7 . 098E+01 

1.  2  64E+02 

6.  444E+02 

1. 3  53E+03 

1200.0 

1 . 903E+00 

7. 618E+00 

6.  876E+01 

1.  224E+02 

6.243E+02 

1 . 311E+03 

1 500. 0 

1 . 853E+00 

7 . 419E+00 

6.  697E+01 

1.  192E+02 

6. 081C+02 

1 . 277E+03 

2000.0 

1 . 824E+00 

7. 301E+00 

6.  590E+01 

1.  173E+02 

5. 984E+02 

1.256E+03 

2  500.  0 

1.823E+00 

7. 297E+00 

6.  586E+01 

1.  172E+02 

5.979E+02 

1. 255E+03 

3000.0 

1 . 833E+00 

7. 339E+00 

6.  623E+01 

1.  179E+02 

6.013E+02 

1. 262E+03 

4000.0 

1 . 867E+00 

7. 472E+00 

6.  743E+01 

1.  200E  +  02 

6. 120E+02 

1 . 285E+03 

5000. 0 

1 . 904E+0  0 

7. 620E+00 

6.877E+01 

1.  224E+02 

6. 240E+02 

1.  310E+03 

6000.0 

1. 940E+00 

7 . 7  64E+00 

7. 006E+01 

1.  247E+02 

6. 3  56E+02 

1 . 334E+03 

7000.0 

1 . 973E+00 

7 . 898E+00 

7. 127E+01 

1 . 269E+Q2 

6.  465E+02 

1 . 357E+03 

8000.0 

2. 004E+0  0 

8 . 0P2E+00 

7 . 238E+01 

1.  28BE+02 

6.  565E+02 

1 . 378E+03 
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Table  2.4 


Range  in  Air  (gm/cra 


2 


) 


MeV/u 

H 

He 

1.000 

2. 907E-03 

3 . 360E-03 

1.200 

3 . 831 E-0  3 

4. 369E-G3 

1 .  500 

5.  561E-03 

fa. 09  5E-03 

2.000 

8.919E-03 

9 . 493E-03 

2.500 

1. 295E-02 

1 . 3  50E-02 

3.000 

1. 762E-02 

1 .809E-02 

4.000 

2. S90E-02 

2 . 896E-02 

5.000 

4.  2  54E-02 

4. 188E-Q2 

6.000 

5.  838E-02 

5. 67  5E-02 

7.000 

7. 640 E -02 

7 . 3  57E-02 

a. ooo 

9. 675E-02 

9. 233E-02 

10.000 

1 . 437E-01 

1. 3  55E-01 

12.000 

1 . 988E-01 

1 . 864E-01 

15.000 

2.961E-01 

2.771E-01 

20.000 

4. 953E-01 

4. 677E-01 

25.000 

7 . 382E-01 

7 . 063E-01 

30  000 

1. 023E+00 

9 . 889E-01 

40.000 

1.716E+00 

1 . 676E+00 

50.000 

2.  562E+0O 

2. 516E+00 

60.000 

3.  554E+00 

3. 500E+00 

70.000 

4. 68 5E+00 

4. 622E+00 

80.000 

5. 948E+00 

5.875E+00 

100.000 

8. 847E+00 

8 . 7  52E+00 

120.000 

1. 221E+01 

1.209E+01 

150.000 

1.S06E+01 

1.7Q9E+01 

200.000 

2.969E+01 

2.943E+01 

250.000 

4. 336E+01 

4. 299E+01 

300.000 

5.877E+01 

5.828E+01 

400.000 

9. 386E+01 

9. 310E+01 

500.000 

1. 335E+02 

1 . 324E+02 

600.000 

1. 765E+02 

1.751E+02 

700.000 

2. 223E+02 

2. 205E+02 

800.000 

2. 701E+02 

2. 679E+02 

1000.000 

3. 703E+02 

3. 673E+02 

1200.000 

4. 747E+02 

4.  709E+02 

1500.000 

6.  360E+0  2 

6.  309E+02 

2000.000 

9. 107E+02 

9 . 034E+02 

2500.000 

1. 187E+03 

1.  178E+03 

3000.000 

1 . 463E+0  3 

1 . 451E+03 

4000. 000 

2. 008E+0  3 

1 . 992E+03 

5000.000 

2.  543E+03 

2.  523E+03 

6000.000 

3. 067E+0  3 

3.043E+03 

7000.000 

3.  583E+03 

3. 554E+03 

8000.000 

4. 009E+03 

4. 056E+03 

C 

0 

Ar 

2. 130E-03 

1 . 972E-Q1 

2. 124E-03 

2. 581E-03 

2. 360E-03 

2. 462E-03 

3. 313E-03 

2. 961E-03 

2 . 981E-03 

4. 684E-03 

4. 124E-03 

3 . 882E-03 

6. 241E-03 

5. 4Q1E-03 

4.837E-03 

7 . 982E-03 

6.813E-03 

5. 849E-Q3 

1 . 202E-02 

1 . OC4E-02 

8. 049E-O3 

1. 677E-02 

1. 378E-02 

1 .048E-02 

2. 220E-02 

1 . 802E-02 

1 .  31  5E-02 

2.83  5E-02 

2. 277E-02 

1 .  60  5F.-02 

3. 518E-02 

2.802E-02 

1 .918E-02 

5. 091E-02 

4. 002E-02 

2 . 612E-02 

6.932E-02 

5. 398E-02 

3.39  5E-02 

1.01QE-01 

7 . 8  54E-02 

4.73  6E-02 

1. 683E-01 

1.287E-01 

7 . 4O2E-02 

2 . 492E-01 

1 . 897E-01 

1.059E-01 

3 . 439E-01 

2. 610E-01 

1.428E-01 

5.  731E-01 

4. 331E-01 

2. 307E-01 

8. 532E-01 

6. 432E-01 

3. 364E-01 

1 .181E+00 

8.893E-01 

4. 594E-01 

1. 555E+00 

1 .  170E+00 

5. 989E-01 

1 , 973E+00 

1 . 482E+00 

7 . 542E-01 

2 . 931E+00 

2. 200E+00 

1.110E+00 

4. 043E+00 

3. 033E+00 

1. 521E+00 

5 . 974E+00 

4. 479E+00 

2 . 234E+00 

9.814E+00 

7 . 3  55E+00 

3.650E+00 

1 . 433E+01 

1. 073E+01 

5.  311E+00 

1 .941E+01 

1. 454E+01 

7 . 180E+00 

3 . 099E+01 

2. 321E+01 

1.143E+01 

4. 406E+01 

3. 299E+01 

1. 623E+01 

5. 827E+01 

4.  362E+01 

2.144E+01 

7.33  5E+01 

5.490E+01 

2. 697E+01 

8.912E+01 

6. 670E+01 

3.275E+01 

1 . 22  2E+02 

9. 143E+01 

4. 486E+01 

1.  566E+02 

1.  172E+02 

5.747E+01 

2.098E+02 

1.  570E+02 

7.69  5E+01 

3. 004E+02 

2.  248E+02 

1 . 1 01E+02 

3.916E+02 

2. 930E+02 

1.436E+02 

4. 826E+02 

3.  611E+02 

1.769E+02 

6. 624E+02 

4. 9  56E+02 

2.428E+02 

8. 388E+02 

6.  276E+02 

3 . 074E+02 

1 . 012E+03 

7. 571E+02 

3 . 708E+02 

1. 182E+03 

8. 843E+02 

4. 331E+02 

1. 349E+03 

1. 009E+03 

4 . 944E+02 

Fe 

2. 204E-03 

2.  52  5E-03 
3. 007E-03 
3 . 824E-03 
4. 666E-03 
5.  540E-03 
7. 393E-03 
9. 392E-03 

1.  154E-02 
1 . 384E-02 
1 . 628E-02 
2. 162E-02 
2. 755E-02 
3.753E-02 
5 . 700E-02 
7 . 991E-02 
1. 061E-01 
1. 679E-01 
2. 413E-01 
3.256E-01 
4. 208E-01 
5. 263E-01 
7. 667E-01 
1 . 044E+00 
1  .  523E+00 
2. 472E+00 

3.  583E+00 
4.832E+00 
7.  67  2E+00 
1 . 087E+01 
1 . 434E+01 
1 . 803E+01 

2.  188E+01 
2 . 994E+01 
3.834E+01 
5. 131E+01 
7 . 341E+01 
9.  566E+01 
1. 179E+02 
1 . 617E+02 
2 . Q48E+02 
2 . 470E+02 
2 . 88 5E+02 
3 . 293E+02 
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2.5  Flax  at  the  Top  of  the  Atmosphere 

The  cosmic  ray  flux  at  the  top  of  the  atmosphere  has  been  measured 
extensively.  A  compilation  of  these  results  and  procedures  for 
estimating  the  fluxes  are  given  in  Adams  et  al.  (1981).  These  fluxes  are 
dependent  upon  various  solar  system  and  near-earth  environmental 
factors.  Among  these  the  11-year  solar  cycle  causes  variations  of  up  to 
an  order  of  magnitude  in  fluxes.  They  are  greatest  at  solar  minimum,  the 
most  recent  of  which  was  late  In  1977.  Solar  flares,  which  are  of  short 
duration  and  occur  at  unpredictable  intervals,  can  cause  increases  of 
many  orders  of  magnitude  in  the  low  energy  flux.  The  magnetic  field  of 
the  earth  however  can  deflect  all  low  energy  charged  particles  below  some 
cutoff  which  can  be  as  high  as  8  GeV/nucleon  (A/Z  =  2)  and  averages  1  or 
2  GeV/nucleon. 

The  cosmic  ray  fluxes  have  errors  of  30  percent  or  more.  This  error 
contributes  to  the  overall  estimate  of  upset  probabilities,  and  is  the 
most  important  error  at  high  altitudes. 

Isotopic  abundances  are  not  critical  in  estimating  soft  upset  rates. 
They  were  taken  from  current  data  as  mentioned  in  Section  2.2. 

2.  6  Results  of  Propagation 

Vertical  flux  of  several  groups  of  cosmic  rays  as  computed  by  the 
methods  of  this  section  are  shown  in  Figure  2.2.  They  are  compared  with 
data  of  Webber  et  al .  (1967).  Good  agreement  with  these  data  is  found. 


3.0  Cosmic  Ray  Fluxes  in  the  Atmosphere 
3.1  Integrated  Flux  in  the  Atmosphere 

In  this  section  the  combination  of  differential  energy  fluxes  into  a 
total  particle  flux  at  any  point  in  the  atmosphere  is  discussed.  The 

particle  flux  at  any  point  is  both  zenith  and  azimuthal  angle  dependent. 

The  two  factors  which  determine  the  flux  in  any  direction  are  geomagnetic 
cutoff  and  pathlength.  The  zenith  angle,  e,  is  the  angle  the  cosmic  ray 
trajectory  makes  with  the  vertical.  Pathlength  in  the  atmosphere  depends 
only  on  this  angle.  We  may  define  a  function  X(o,x0)  (see  Section  3.2) 

giving  the  number  of  grams  to  the  top  of  the  atmosphere  proceeding  at  a 

zenith  angle  <5  from  a  point  whose  vertical  depth  in  the  atmosphere  is 
x0  g/cm2.  This  function  has  the  property: 

X(0,Xq)  =  x0  (3.1-1) 

The  geomagnetic  cutoff  at  a  point  depends  both  on  the  zenith  and 
azimuthal  angles.  We  define  the  azimuthal  angle,  <t> ,  as  the  incoming 
direction  of  the  cosmic  ray  as  measured  with  respect  to  geomagnetic 
east.  At  any  geographic  point,  the  cutoff  vs.  zenith  and  azimuthal 
angles  is  completely  determined  by  the  vertical  cutoff.  Thus,  the  lower 
limit  of  rigidities  allowed  within  the  atmosphere  at  any  place  is  given 
by  a  function  R(e,0,Rv)  (see  Section  3.3).  The  rigidity  is  related  to 
energy/nucleon  for  a  particular  species  i  by: 
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R(GV/c) 


(3.1-2) 


A1  ^E(E  +  1863) 

ioooz1 


where  E  is  expressed  in  MeV/nucleon. 

Given  the  pathlength  and  rigidity  cutoff  functions  and  the  flux  vs. 
energy  and  pathlength,  J(x,E),  as  computed  in  Section  2,  the  total  flux 
at  a  point  is: 


F  (x  ,E) 
T  0 


7T 

j  sin©  de  J' 


2n 

d 0  J(X(e,xo),E)  g(e,0,E) 


(3.1-3) 


where  J( X , E )  is  the  calculated  particle  flux  (see  Section  3)  for  energy  E 
and  pathlength  X. 

The  function  g(e,0,E)  *  1  if  the  particle  was  able  to  enter  the  top 
of  the  atmosphere  and  =>  0  otherwise.  We  evaluate  g(«,0,E)  by  comparing 
the  range  of  the  nuclide  at  energy  E  in  air  with  the  range  at  cutoff 
rigidity  minus  the  pathlength.  If 

Range(E)  ^  Range(F.cutoff)  -  X(e,x0)  (3.1-4) 

then  the  range  at  the  top  of  the  atmosphere  was  greater  than  the  range  at 
cutoff  energy  (Implying  the  energy  was  greater  than  cutoff  energy),  and 
g(©,0,E)  -  1. 


The  flux  calculations  are  performed  for  a  limited  number  of 
pathlengths  (In  this  work,  7).  The  flux  at  any  other  pathlength  is 
determined  by  Interpolating  between  these  fluxes  (or  extrapolating  if 
X(e,x0)  >  100  g/cm2.  The  fluxes  at  two  nearby  points  x0  and 

are  J0  and  Jl  respectively.  Assuming  an  exponential  decrease  in  flux 
with  pathlength  the  flux  at  x  g/cm2  is: 


J(x) 


JQ  exp 


j  1 0(0^) 
(  i*i-  V 


(3.1-5) 


Figure  3.1  shows  total  flux  levels  for  the  altitudes  of  interest 
in  this  report  versus  geomagnetic  cutoff.  From  the  figure  one  sees 
that  total  flux  can  be  reduced  more  than  a  factor  ten  by  the 
geomagnetic  cutoff.  Roughly  equal  spacing  of  the  curves  corresponds 
to  the  roughly  exponential  drop  in  flux  with  pathlength. 


3.2  Zenith  Angle  Dependence  of  Atmospheric  Depth 

The  variation  of  pathlength  with  zenith  angle  and  altitude  is 
computed  in  this  section.  Figure  3.2  shows  a  schematic  diagram  of 
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earth's  atmosphere  at  some  altitude  h.  The  vertical  distance  to  the 
top  of  the  atmosphere  is  denoted  by  x.  We  wish  to  calculate  the 
distance,  y,  to  the  top  of  the  atmosphere  at  a  zenith  angle  9. 
Applying  the  law  of  cosines  we  have: 

y2  +  2y(Re+h)  cos  e  +  (Re+h)2  =  (Re+h  +  x)2  (3.2-1) 

therefore: 


(Re+h) 


1/2 


COS  9 


The  square  root  may  be  re-expressed  as: 


cos  e 


n  + — £l 

L  (Rfl+h)c 


h)cos^o 


/  2 


(Re+h')2cos23l 


(3.2-2) 


(3.2-3) 


In  this  work,  the  altitudes  under  consideration  are  less  than  64  km 
whereas  the  radius  of  Earth  Is  -  6400  km.  The  third  term  In  brackets 
may  be  neglected,  while  so  long  as  9  <  80°  the  remainder  can  be 
expanded  to  yield 


y  =  x  sec  9 


(3.2-4) 


Therefore,  for  zenith  angles  which  encompass  the  majority  of  cosmic  rays 
at  a  point,  the  relation  between  y  and  x  is  linear. 

The  linear  relationship  between  x  and  y  implies  a  linear 
relationship  between  the  pathlength  along  the  vertical  to  the  pathlength 
at  9.  In  moving  vertically  between  x  and  x  +  dx,  a  density  0(x)  is 
encountered  and  a  pathlength  p(x)dx  is  traversed.  The  corresponding 
move  at  an  angle  9  Is  over  a  distance  dy  =  dx  sec  9  and  therefore 
Involves  traversal  of  precisely  sec  9  times  the  vertical  pathlength  - 
Independent  of  x.  Tnus,  the  pathlength  vs.  zenith  angle  relation  Is 
simply 


X(o,x0)  =  x0  sec  9  (3.2-5) 

At  the  point  where  this  relation  breaks  down  (e  =»  80°)  the  cosmic  ray 
flux  is  negligible,  having  passed  through  almost  6  times  as  much  air  as 
the  vertical  flux. 


3.3  Geomagnetic  Cutoff 

Earth's  magnetic  field  serves  as  an  extremely  effective  shield  of 
low  to  medium  energy  cosmic  rays.  This  shielding  action  is 
characterized  by  a  vertical  rigidity  cutoff  at  each  geographical 
location  (see  energy/ rigidity  relation,  eq.  3.1-2).  Positively  charged 
particles  having  rigidities  below  the  cutoff  cannot  arrive  at  a  point 
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vertically  or  from  the  east  because  they  are  deflected  away.  Arriving 
from  the  west  they  are  deflected  toward  Earth  and  hence  are  able  to 
penetrate  at  lower  energies. 

Vertical  geomagnetic  cutoffs  have  been  computed  at  20  km  altitude  by 
Shea  and  Smart  (1975).  For  roughly  800  geographical  positions,  charged 
particles  arriving  vertically  were  propagated  backward  through  a  precise 
model  of  the  geomagnetic  field.  Such  trajectories,  when  the  rigidity  Is 
below  cutoff,  eventually  Intersect  the  earth's  surface.  When  the 
rigidity  Is  above  cutoff  they  escape  earth.  This  method  gives  an 
accurate  determination  of  the  cutoff  rigidity.  A  contour  plot  of  the 
vertical  cutoff  data  Is  presented  in  Figure  3.3.  Note  that  more  than 
half  earth's  surface  has  a  cutoff  rigidity  of  greater  than  4  GV.  Since 
the  average  cosmic  ray  energy  is  roughly  1  or  2  GeV/nucleon  most  cosmic 
rays  are  deflected  from  earth  by  its  magnetic  field, 

The  cutoff  rigidity  in  directions  other  than  the  vertical  depends  on 
the  angle  of  arrival  (see  Fig.  3.4).  Low  energy  particles  arrive 
predominantly  west,  while  particles  with  energy  much  greater  than  the 
cutoff  arrive  from  any  direction.  Exact  arrival  angles  are  described  by 
the  Stormer  cones.  The  relation  between  the  5tormer  cone  angle,  y,  and 
cutoff  is  approximately 

o  _ ^ _ __  (3.3-1) 

[1  +  (1  -  cos  V'cos^  x)  ^ 

where 


cos4  \  *  Rv  (l  +  altitude  (km) 

T775  - 5371 — 

Plots  of  the  ratio  of  cutoff  to  vertical  cutoff  for  various  vertical 
cutoffs  and  Stormer  cone  angles  are  displayed  in  Fig.  3.4. 

The  Stormer  cone  angle,  y»  is  related  to  the  azimuthal  and  zenith 
angles  by 


): 


(3.3-2) 


cos  y  -  cos  (6  sin  e  (3.3-3) 

Eq.  3.3-1  through  3.3-3  allow  the  geomagnetic  cutoff  for  any  arriving 
angle  to  be  determined  from  the  vertical  cutoff  and  altitude.  Further 
discussion  of  the  geomagnetic  field  effect  on  cosmic  rays  is  contained 
in  Adams  et  al.  (1983). 


4.0  Single  Event  Upsets 
4.  1  LET  Spectra 

Up  to  this  point  we  have  discussed  the  calculation  of  cosmic  ray 
spectra  in  the  atmosphere  for  each  charge  species.  There  is  no  direct 
relation  between  these  spectra  and  single  event  upsets  because  the 
rate  of  energy  deposit  by  a  cosmic  ray  bears  a  complicated  functional 
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relation  to  its  total  energy.  The  LET  (linear  energy  transfer)  or 
stopping  power  of  a  charged  particle  is  its  rate  of  energy  loss  in 

matter.  Table  4.1  shows  the  stopping  power  in  silicon  for  several 
elements  in  MeV/(g/cm2}.  It  is  notable  that  the  stopping  power  is 

at  a  minimum  and  roughly  constant  above  1  SeV/nucleon. 

The  transformation  from  energy  spectrum  to  LET  spectrum,  is: 

J(L)  -  E01  (E)  (4.1-1) 

where  the  sum  over  species  is  necessary  because  the  stopping  power  Is 
charge  dependent.  The  differential  LET  spectrum  represents  the  energy 
transfer  characteristics  of  all  species  in  one  spectrum.  Because  L  is 
not  a  monotonic  function  of  E  Its  derivative  can  be  zero,  thus  J(L) 
has  singular  points  corresponding  to  stopping  power  minima  and  maxima 
for  all  the  elements.  These  singularities  are  avoided  in  numerical 
applications  by  using  LET  and  energy  bins,  and  defining  the  LET 

spectrum  by 

J(L)  -EJ1  (E)  $  (4.1-2) 

In  our  work  we  have  divided  the  LET  Interval  (1  <_  L  <  105 

MeV/(g/cm2))  and  the  energy  Interval  [1  <  E  <  8  GeV/nucleonf  into 

500  bins  with  special  treatment  of  partlcler above  8  GeV/nucleon. 

In  practice  It  is  convenient  to  use  the  integral  LET  spectrum 

oo 

N(L)  -  f  J(L)  dL  (4.1-3) 

.  L 

because  it  does  not  fluctuate,  but  decreases  monotonical ly  from 
L  «  0.  Integral  LET  spectra  for  cosmic  rays  at  various  altitudes  and 
geomagnetic  cutoffs  are  shown  in  Fig.  4. 1-4, 9.  Interpolation  allows 
the  spectrum  to  be  found  at  any  altitude  between  50  and  150  thousand 

feet  at  any  location.  For  example,  Fig.  4.10  shows  the  integral  LET 

spectrum  at  103  MeV/(g/cm2)  for  all  altitudes  and  cutoffs. 

Figure  4.1  shows  the  typical  behavior  of  the  LET  spectrum  as  a 
function  of  altitude.  At  150,000  ft.,  where  very  little  slowing  down 
has  taken  place,  the  geomagnetic  field  essentially  eliminates  Ions 
with  LET  >  2  x  103  MeV/(g/cm2).  This  occurs  because  all  heavily 

ionizing  (i.e.,  low  energy)  Fe  are  deflected  while  higher  charge 
species  are  extremely  rare.  By  85,000  ft.  many  high  energy  Fe  ions 
have  slowed  down  because  of  Ionization  loss  and  the  spectrum  is 

appreciable  to  10*4  MeV/(g/cm2). 

4.2  Soft  Upsets 

The  relation  between  soft  upsets  and  the  LET  spectrum  is  governed 
by  the  critical  charge  of  the  device,  Qc,  and  the  geometry  of  the 
sensitive  volume.  An  energy  deposition  in  the  sensitive  volume  causes 
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Table  4.1 


Stopping  Power 


in  Silicon  (MeV/ (grn/cra 


2 


)> 


MeV/u 

H 

He 

1.000 

1.7  59E  02 

7 . 046E 

02 

1.200 

1.  57 6E  02 

6.305E 

02 

1. 500 

1.369E  02 

5.466E 

02 

2.000 

1.133E  02 

4. 502E 

02 

2.  500 

9.710E  01 

3.846E 

02 

3.000 

8. 524E  01 

3. 366E 

02 

4.000 

6. 873E  01 

2.704E 

02 

5.0C0 

5.  892E  01 

2. 308E 

02 

6.000 

5.  1 31E  01 

2.006E 

02 

7.000 

4.  5  58E  01 

1.780E 

02 

8.000 

4.  HOE  01 

1. 604E 

02 

10.000 

3.452E  01 

1. 347E 

02 

12.000 

2.992E  01 

1.170E 

02 

15.000 

2.  51 4E  01 

9.888E 

01 

20.000 

2.014E  01 

a. 006E 

01 

25.000 

1.69  5E  01 

6.775E 

01 

30.000 

1.469E  01 

5.880E 

01 

40.000 

1.172E  01 

4.693E 

01 

50.000 

9.854E  00 

3.944E 

01 

60.000 

8. 561E  00 

3.426E 

01 

70.000 

7. 612E  00 

3.046E 

01 

80.000 

6.883E  00 

2.755E 

01 

100.000 

5.837E  00 

2. 336E 

01 

120.000 

5.  120E  00 

2.049E 

01 

150.000 

4.  38 5E  00 

1.7  55E 

01 

200.000 

3.631E  00 

1.453E 

01 

250. 000 

3.170E  00 

1. 269E 

01 

300.000 

2.860E  00 

1. 145E 

01 

400.000 

2.472E  00 

9 . 89  5E 

00 

500. 000 

2.242E  00 

8.976E 

00 

600.000 

2.094E  00 

8.380E 

00 

700.000 

1.991E  00 

7.971E 

00 

800.000 

1.918E  00 

7. 677E 

00 

1000.000 

1.818E  00 

7.277E 

00 

1200.000 

1.758E  00 

7.038E 

00 

1500.000 

1.708E  00 

6.837E 

00 

2000.000 

1. 67 4E  00 

6. 702E 

00 

2500.000 

1 . 667E  00 

6. 673E 

00 

3000. 000 

1.671E  00 

6. 689E 

00 

4000.000 

1. 69 IE  00 

6.770E 

00 

5000. 000 

1.716E  00 

6. 863E 

00 

6000.000 

1.740E  00 

6.  965E 

00 

7000.000 

1.763E  00 

7. 0  55E 

00 

8000. 000 

1.784E  00 

7. 139E 

00 

C  0  Ar 

4.401E  03  6.768E  03  1 . 890E  04 

4.1CBE  03  6.492E  03  1 • 894E  04 

3.851E  03  6.089E  03  1.865E  04 

3.409E  03  5.493E  03  1 . 790E  04 

3.056E  03  4.993E  03  L.706E  04 

2.769E  03  4.  575E  03  1.623E  04 

2.330E  03  3.916E  03  1.472E  04 

2.032E  03  3.453E  03  1 .  348E  04 

1.794E  03  3.078E  03  1. .  242E  04 

1 . 607E  03  2. 779E  03  1. 151E  04 

1.457E  03  2.  534E  03  1.074E  04 

1.232E  03  2.158E  03  9.473E  03 

1.070E  03  1.0G3E  03  8.487E  03 

9.000E  02  1.588E  03  7.36XE  03 

7.209E  02  1.274E  03  6.051E  03 

6.082E  02  1.076E  03  5.169E  03 

5.2S7E  02  9.370E  02  4. 542E  03 

4.225E  02  7.502E  02  3.703E  03 

3.552E  02  6.310E  02  3.146E  03 

3.087E  02  5.486E  02  2 . 748E  03 

2.745E  02  4.879E  02  2.452E  03 

2 . 483 E  02  4.414E  02  2.223E  03 

2.106E  02  3.745E  02  1.891E  03 

1.847E  02  3.286E  02  1.663E  03 

1.583E  02  2.816E  02  1.427E  03 

1 . 3.1  IE  02  2.  333E  02  1.1B5E  03 

1.145E  02  2.037E  02  1.036E  03 

1.0332  02  1.838F,  02  9.355E  02 

8 . 93 0E  01  1.590E  02  8.097E  02 

8.102E  01  1.442E  02  7.351E  02 

7.56SE  01  1.347E  02  6.867E  02 

7.195E  01  1.281E  02  6.  533E  02 

6 . 93  IE  01  1.234E  02  6.294E  02 

6.S69E  01  1.170E  02  5.968E  02 

6.354E  01  1.131E  02  5.773E  02 

6.173E  01  1.099E  02  5.609E  02 

6.0S1E  01  1.077E  02  5.498E  02 

6.025E  01  1 . 073 E  02  5.474E  02 

6.039E  01  1.075E  02  5.487E  02 

6 , 11 IE  01  1.088E  02  5. 551E  02 

6.  199E  01  1.104E  02  5. 630E  02 

6.286E  01  1.119E  02  5.709E  02 

6.368E  01  1.134E  02  5.782E  02 

6.444E  01  1.147E  02  5.350E  02 


Fe 

2  .  7  7  3E 

04 

2.816E 

04 

2.835E 

04 

2.802E 

04 

2 . 733E 

04 

2. 650E 

04 

2. 477E 

04 

2. 319E 

04 

2 . 174E 

04 

2. 046E 

04 

1 .931E 

04 

1.737E 

04 

1. 580E 

04 

1. 393E 

04 

1 . 170E 

04 

1 . 012E 

04 

8. 969E 

03 

7 . 408E 

03 

6.  377E 

03 

5.  61 GE 

03 

5. 03QE 

03 

4. 583E 

03 

3.913E 

03 

3.455E 

03 

2.974E 

03 

2.47  5E 

03 

2. 167E 

03 

1.959E 

03 

1. 697E 

03 

1.  542E 

03 

1. 441E 

03 

1 . 371E 

03 

1. 32 2E 

03 

1. 2  53E 

03 

1. 213E 

03 

1. 178E 

03 

1.  155E 

03 

1. 150E 

03 

1. 153E 

03 

1. 166E 

03 

1 . 182E 

03 

1. 199E 

03 

1. 214E 

03 

1. 228E 

03 
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free  electron-hole  pairs  to  be  created  at  a  rate  of  1  pair  per  3.6  eV 
deposition.  This  charge  is  then,  by  definition  of  the  sensitive 
volume,  able  to  accumulate  and  cause  a  change  in  device  state.  Device 
sensitivities  are  discussed  in  a  number  of  articles  (Pickel  and 
Blandford,  1980;  Petersen  et  al.,  1982;  Ziegler  and  Lanford,  1979). 

The  critical  charge  is  a  few  pC  for  devices  with  length  scales  on  the 
order  of  10-20  ym.  It  varies  roughly  as  the  inverse  square  of  the 
length  scale.  Radiation  hardened  devices  have  a  critical  charge  about 
twice  as  large.  The  total  charge  freed  when  & E  is  deposited  is 

q(pO  =  A-T?:r}  (4,1“4) 

Thus  the  upset  rate  Is  the  number  of  particles  per  unit  time  which 
deposit  an  energy  greater  than  or  equal  to 

*EC  -  22.5  Qc  (4.1-5) 

Since  ions  lose  energy  linearly  at  a  rate  equal  to  their  LET  (in  thin 
slabs  of  matter),  the  upset  rate  depends  only  on  the  LET  spectrum  and  the 
distribution  of  chord  lengths  in  the  sensitive  volume. 


The  chord  length  distribution  C(p)  Is  a  probability  density. 
Assuming  the  sensitive  volume  is  subject  to  an  isotropic  cosmic  ray  flux, 
C(p)  dp  is  the  probability  that  a  cosmic  ray  will  pass  through  between  p 
and  p  +  dp  of  the  region.  This  is  a  purely  geometric  quantity  which  is 
known  exactly  for  any  rectangular  parallelepiped  (Pickel  and  Blandford, 
1980).  The  energy  deposited  by  a  charged  particle  over  the  length  p  is 

aE  -  pLp  (4.1-6) 


where  p  is  the  density  of  silicon  (2.32  g/cm3). 
any  charged  particle  with  LET  greater  than 

i  _  22,5  q. 


For  a  chord  length  p 
(4.1-7) 


causes  a  device  upset.  The  total  number  of  upsets/sec  is  therefore 

00 

U  “  \  /N(Lm1n)  c(p)  dp,  (4-1_8) 

o 

where  Ap  is  the  mean  projected  area  of  the  sensitive  volume,  equal  to 
one  quarter  the  surface  area  of  a  rectangular  parallelepiped,  and  C(p)  is 
the  chord  length  distribution.  N(Lm-jn)  is  given  by  Equation  4.1-3. 


Figures  4.11  through  4.14  show  typical  upset  rates  and  their 
variation  with  parameters  of  the  calculation.  In  Figure  4.11  the 
variation  with  altitude  show  that  upset  rates  for  a  sensitive  volume 
5  x  10  x  10  ym  and  critical  charge  between  .1  and  1  pC  are  strongly 
dependent  on  atmospheric  depth.  At  the  top  of  the  atmosphere  the  cutoff 
excludes  highly  ionizing  Fe,  while  at  lower  altitudes  ionization  loss 
restores  the  low  energy  iron  flux.  At  75,000  ft.  the  variation  with 
cutoff  is  shown  In  Figure  4.12.  Highly  Ionizing  Fe  are  restored  from  a 
cutoff  of  1  GV,  but  not  at  the  4  GV  level.  The  breakdown  into  element 
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groups  (Figure  4.13)  shows  that  for  critical  charges  above  0.02  pC,  the 
heaviest  ions  contribute  most  to  soft  upsets.  Lighter  ions  are  of 
Importance  in  proportion  to  their  charge.  Figure  4.14  shows  the 
variation  of  several  orders  of  magnitude  possible  with  varying  device 
sensitivity. 


5.0  Computer  Code 

In  this  section  the  computer  code  developed  to  estimate  soft  upset 
rates  in  the  atmosphere  is  presented.  These  programs  may  be  classified 
as  follows: 

(i)  Data  Filer  vj  Data  File  Preparation  Programs 

( i  1 )  Air  Prop-.  hn  Programs 

(Hi)  Post-proce-.  .ng  Programs 

(1v)  Command  File  Generator 

Classes  (1)  through  (111)  are  VAX-11  FORTRAN  programs  and  data 
files.  This  FORTRAN  version  Is  an  extension  of  the  standard  FORTRAN-77 
and  the  programs  may  not  be  compatible  with  other  systems.  The  structure 
of  this  system  of  programs  and  data  files  is  shown  in  Figure  5.1.  The 
Command  File  Generator  Is  written  in  the  VAX/VMS  command  language,  DCL 
(Version  3.2).  This  file  generates  a  command  file  which  associates  the 
programs  and  data  files  correctly  and  requests  all  required  input  data. 
The  command  file  may  be  submitted  for  execution  when  convenient, 

5.1  Data  Files  and  Data  File  Preparation 

There  are  four  input  data  files  required  for  execution  of  these 
programs. 

1)  SILPOW  -  Stopping  powers  of  H,  He,  C,  0,  Ar,  and  Fe  in  Si  (same 
as  Table  IV. 1).  Format  (F8.3,6(3X,1PE9.3) ) . 

2)  ASTPOW  -  Stopping  powers  of  H,  He,  C,  0,  Ar  and  Fe  in  air  (same 
as  Table  II. 3).  Format  (F8.3,6(3X,1PE9.3) ) . 

3)  ARANGE  -  Range  of  H,  He,  C,  0,  Ar  and  Fe  in  air  (same  as  Table 
II. 4)  Format  (F8.3,6(3X,1PE9.3) ) . 

4)  S0URCE3  -  List  of  Isotopic  Species  (same  as  Table  11.2),  Format 
(5(13, A2, 13) ,F6.3) . 

In  addition,  four  data  file  preparation  programs  create  the  remaining 
four  files: 

5)  QAIR3  -  List  of  partial  cross  sections  for  production  of  lighter 
nuclides  from  heavier  nuclides  in  collision  with  air,  produced  by  the 
program,  QAIR,  which  requires  the  subroutine  YIELDX  (IZ,  IA,  OZ,  JA,  E, 
S).  YIELDX  returns  the  partial  cross  section,  S,  for  the  process  (N^  + 
p  »  Nj  +  p  +  anything)  at  energy,  E  (MeV/nucleon) . 


18 


6)  FLUXES  -  Flaxes  at  the  top  of  the  atmosphere  for  the  Isotopic 
species  In  SQURCE3.  The  program  FLUX  uses  Adams  et  al.  (1981)  model  of 
cosmic  ray  fluxes  as  modified  by  recent  experimental  work. 

7)  RANGIN  -  Produced  by  RANGER  and  contains  data  on  flux  loss  in 
passing  through  .1  gm/cm2  air. 

8)  POWIN  -  Produced  by  POWER  and  contains  stopping  power  data  on  the 
nuclides  in  S0URCE3,  interpolated  from  ASTPOW. 

The  data  file  preparation  programs  are  listed  in  Appendices  1  through 
9.  A  subroutine,  CUBSPL,  which  performs  cubic  spline  Interpolation  is 
required  for  execution  of  RANGER.  This  subroutine  is  shown  in  Appendix  5. 


5.2  Air  Propagation  Programs 

The  techniques  used  in  the  air  propagation  programs  are  discussed  In 
Section  2.  The  procedure  is  separated  into  two  parts.  In  the  first 
(APROP)  propagation  through  all  but  the  last  gram  of  air  is  performed. 
Only  fragmentation  for  the  last  gram  Is  performed.  In  the  second  (GRAM) 
the  last  gram  of  Ionization  loss  is  performed,  meanwhile  extending  the 

minimum  energy  from  30  MeV/u  to  1  MeV/u.  These  programs  are  shown  in 
Appendices  6  and  7. 

The  two  subroutines  of  APROP  (appended  to  the  program)  are  AIRMAT  and 
ILOSS .  AIRMAT  constructs  the  matrix  of  total  and  partial  cross 
sections.  The  partial  cross  sections  are  read  in  from  QAIR3  while  the 
total  cross  sections  are  calculated  using  2.3-1  and  2.3-2.  In  ILOSS  the 
coefficients  In  2.1-7  pertaining  to  ionization  loss  are  calculated  as 

needed. 

5.3  Post-Processing 

The  three  post-processing  programs  exhibited  In  Figure  6.1  are  found 
In  Appendices  8,  9  and  10.  These  programs  take  the  results  of  air 
propagation  and  convert  them  into  the  quantities  of  interest  in  this 
report. 

ATMOS  interpolates  the  seven  flux  files  at  0,  5,  15,  30,  50,  75,  and 
100  gm/cm2  and  Integrates  them  over  zenith  and  azimuthal  angles  to  get 

the  total  integrated  flux.  These  procedures  are  discussed  in  Section  3. 
The  functions  RANGE  and  STPOW  estimate  the  range-energy  and  stopping 
power  relations  in  air  using  the  Bethe-Bloch  theory.  These  are  used  in 
combination  with  CUTFACT  to  evaluate  the  Inequality  3.1-4.  CUTFACT 
determines  the  angular  dependence  of  the  geomagnetic  cutoff.  The 

subroutine  ALTGRM  computes  the  altitude  vs.  grammage  relation  of  Fig. 
2.1.  FLUX  performs  the  exponential  Interpolation  (Eq.  3.1-5). 

LETINTT  performs  the  conversion  from  flux  to  LET  spectrum  for  all 
species.  It  calls  the  cube  spline  subroutine,  CUBSPL.  Methods  are 
discussed  in  Section  5.1 . 
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UPSET  converts  the  LET  spectrum  Into  an  upset  vs.  critical  charge 
profile.  It  depends  on  geometrical  properties  of  the  sensitive  region 
which  is  assumed  to  be  a  rectangular  parallelepiped.  The  function  DIFPLD 
evaluates  the  analytic  form  of  the  differential  pathlength  distribution. 
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PATHLENGTH  (X)  vs.  ALTITUDE  (A3  BETWEEN 
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Figure  2.1  Altitude  vs.  pathlength  in  the  atmosphere  between 

50,000  ft.  and  150,000  ft.  Pathlength  is  proportional 
to  pressure  and  varies  roughly  exponentially  with 
altitude. 


Figure  2,2  Comparison  of  vertical  cosmic  ray  fluxes  in  the 

atmosphere  (Webber  et  at.,  1967)  with  computations  of 
this  work. 


FLUX  (particles/m2sec) 


INTEGRAL  FLUX  vs. 
CUTOFF  FOR  VARIOUS  ALTITUDES 


Figure  3.1  Integral  cosmic  ray  flux  (mostly  protons)  as  a  function 
of  altitude  and  geomagnetic  cutoff.  Equal  spacing  of 
curves  shows  roughly  exponential  decrease  with 
Increasing  depth  in  the  atmosphere. 
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Figure  3.2  Schematic  diagram  of  relation  between  vertical 

pathlength,  x,  to  a  point  in  the  atmosphere  and  the 
pathlength,  y,  at  zenith  angle  e. 


RATIO  OF  CUTOFF  TO  VERTICAL  CUTOFF 


STORMER  CONE  ANGLE  FOR  POSITIVELY  CHARGED  PARTICLES 

(Degress  from  West) 

Figure  3.4  Cutoff  rigidity  vs.  Stormer  cone  angle  for  various 

vertical  cutoffs.  Widest  variations  in  cutoff  occur  in 
regions  of  high  vertical  cutoff. 
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PARTICLES/m2  sec 


LET  (MeV  cm*  /gm) 


Figure  4.1 


LET  spectra  at  1  GV  cutoff  for  various  altitudes. 
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PARTICLES/m2  sec 


I  10  10*  10s  I04  10° 


LET  (MeV  cm2  /gm) 

Figure  4.2  LET  spectra  at  55,000  ft.  for  several  cutoffs. 
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PARTICLES/m2  sec 


PARTICLES/m2  sec 


Figure  4.4  LET  spectra  at  6b, 000  ft.  for  several  cutoffs. 
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PARTICLES/m2  sec 


Figure  4.5  LET  spectra  at  75,000  ft.  for  several  cutoffs. 
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PARTICLES/m2  sec 


10  10*  10s  I04  10° 

LET  (MeV  cm  2  /gm) 


Figure  4.6  LET  spectra  at  85,000  ft.  for  several  cutoffs. 
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PARTICLES/m2  sec 


LET  <MeV  cm2  /gm) 

Figure  4.7  LET  spectra  at  100,000  ft.  for  several  cutoffs. 
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PARTICLES 


F 


Figure  4.8  LET  spectra  at  125,000  ft.  for  several  cutoffs. 


i 
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Figure  4.9  LET  spectra  at  150,000  ft,  for  several  cutoffs. 


(particles/mz  sec) 


Qc ( pC) 


Figure  4.11  Variation  of  soft  upset  rate  with  altitude. 
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Qc  ( PC) 

Figure  4.12  Variation  of  soft  upset  rate  with  geomagnetic  cutoff. 
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Figure  4.14  Variation  of  soft  upset  rate  with  sensitive  volume. 
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POST  PROCESSING  PROPAGATION  DATA  FILES  AND  DATA  FILE  PREPARATION 


Figure  5.1  Flow  chart  of  routines  and  data  files  used  in  making  the 
calculation  of  this  report. 
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APPENDIX  1 


PRQORAM  0A1R 

4*4#4#4*44#4#*44444#44##444#*#**4##*#*4*4*#4*4#*#**##44444*#44# 

***  THIS  PROORAM  COMPUTES  PARTIAL  CROSS  SECTION  FOR  CREATION 
*##  OF  ONE  ISOTOPIC  SPECIES  FROM  ANOTHER  IN  COLLISIONS  WITH 
#**  AIR  NUCLEI .  IT  IS  BASED  ON  PROTON-NUCLEUS  PARTIAL 
#**  CROSS  SECi-ICHil  (YIELDX.  FOR  >  COMPUTED  FROM  THE  SEMI- 
44*  EMPIRICAL  FORMULAS.  THE  LIST  DF  SPECIES  (B0URCE3.DAT) 

#**  IS  REQUIRED  INPUT.  QAJR3,  DAT  IS  PRODUCED  AS  OUTPUT. 

#4**4#  444*#4#*#444,,444  4#4#4###4#*#*##4##***##*#4#*#4#*4*##***#* 

PARAMETER  N-Vi,  N 

DIMENSION  12(5,  N).  ‘iX(3,  N),  I Y (  3,  N > ,  FI  ( N ) ,  Q  ( N.  N ) ,  81  ( N ) ,  SS  ( N ) 
DIMENSION  M1(N)  ,  H2lN>,  IQfN,  N > ,  QFCN),FZ<26) 

DATA  KZ/7.  E4,  4.  E3,  19.  7,  9.  7,  28.  4,  100.  ■  89.  2,  98.  1,  1.6.  1.4 

1.  2.8,19.0,  3.3.14.0,  0,  9,  3.0,  1.7,  3.1,  2.2,  4.7 

1,  1.  a,  3.  3,  1.  7,  3.  3,  2.  4,  20.  5/ 

DATA  SI,  S2/N*0,  ,  N#0,  / 

OPENtUNIT-1,  FILE"  'SOURCES.  DAT',  BTATUB-  'OLD ' ) 

0PEN(UNIT-2,  FILE-'QAIR3.  DAT  ',  STATUS-  'NEW ' ) 

C  #**4  4*#*4##**#4######*##4######4**##44#'*44  #*##**#*4*##***####*4 

C  #**  READ  IN  LIST  OF  SPECIEB. 

C  4*444#4##4*4*###444###4#44###******4#**##*###*#*44#*##4##*#**4* 

READ  (1,1)  (  <  IZ  <  J,  I),  IX(J,  1),  IY(J,  I>,  J-l,  9),QF(I),  I-1,N) 

1  FORMAT  ((5(13,  A2,  13), F6.  3)) 

C  #*#44*#*##**#*4*4###4#4*#***44*****4#*#-»4*##4*  **4*4#**  *'##****#* 

C  *#*  FOR  EACH  TARSET  SPECIEB  I,  COMPUTE  THE  CR08B  SECTION  TO 

C  ##*  CREATE  THE  PRODUCT  SPECIES  J  AND  UP  TO  4  OTHER  NUCLIDES 

C  ***  WHICH  ARE  ORQUPSD  WITH  J  IN  THE  SPECIES  LIST. 

C  444444  444444  4 444 44 #4 44 4 4 44 44 44 4444 4444 44 44 4444 4444 44# 444 44 44 444 

DO  I-l.NI 
II-IZU,  I) 

Fill  )-FZ(II)»QF(I) 

DO  J-l,  NI 
Q<  J,  I )  -0. 

DO  K-l,  S 

IF  <IY<K,  Jl.Efl.  0)  00  TO  18 

Z1-IZ(  1,1) 

At-IY(  1,1) 

Z  -IZ(K.J) 

A  -IY<  K,  U> 

CALL  YIELDX  (IZU,  I),  I  Y(l,  I),  IZ(K,  J).  lYtK,  J),  2300.  ,  QJ) 

ED-1. 

IF  ( IZ  <K,  J).  OT,  3  .AND.  A,  LT.  Al/2.  )  ED-3.  *EXP  ( -8.  #A/Ai  ) 

EL-1.  +  4*<  1.  +.  02* (  Zl/Z>#*8)#(  l.  -1.  94Z/Z1  ) 

IF  <  IZ  <K,  J).  LT,  3  .OR.  I  Z  !K,  J ) ,  OT.  5  )  EL-1. 

Q(  J,  I )  -0(  J,  I)+QJ#2.  #ED»EL 
END  DO 

12  SI  (I)-Sl  ( I)#a<  J,  I ) 

32  ( J)  -S2  <  U)  +G<  J,  I ) 

END  DO 

C  #44444  444444  4444  444444*4  4 4*4444*4* 4*4# 444 44444#,,  44444444*4*4  44* 

C  4*4-  PUT  IN  INTEOER  FORM. 

0  44 4# 44  444444  44 4  4 44 #4  4 444  44 44 44 4 *44 4 4 44 4  44444 4 4 44 4 4  44 4**4  444  ^#44 

DO  J-l.  NI 

IO(J,  I)-U(U,  I!+.  5 
END  DO 

C  44 4444 444444 4444 44444444 44 44 4*4444 44444 4444*4 444 #444 44*4 4444*44 

C  ###  OUTPUT  TO  QAIR3.  DAT 

C  4*444* *4**44 *44* #4*44*4* 44*4*4*4*4*4****** #444 4444444*****4*4** 

WRITE(2,  2)  (IZ  ( K,  I),  IX(K,  I),  IY(K,  I ) ,  K-l ,  5 ) ,  F I  (  I  ) , 

S<  (IQ(  J,  I).  J-l.NI) 

2  FORMAT  <S<  13,  A2,  IS),  30X,  FIO.  2/4(2014/),  1614) 

END  DO 

END 
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program  flux 

####■#*  4######*  *■##*•*  4#4#*#*###4#####*#»*4###**#**#***#*#*#*##**4 

*##  THIS  PROGRAM  CREATES  AN  INPUT  FILE  -  FLUXES,  DAT  -  FOR  USE 
4##  IN  THE  AIR  PROPAGATION  PROGRAM.  THE  FUNCTION  CRFE 
*#*  CONTAINS  THE  COSMIC  RAV  FLUX  MODEL  FROM  NRL  MEMORANDUM 
**#  REPORT  4306  i COSMIC  RAY  EFFECTS  DM  MICROELECTRONICS, 

•***  PART  I).  THIS  PROGRAM  ALSO  CREATES  A  FILE  -  QUTFILEO.DAT 
#4#  -  WHICH  CONTAINS  THE  INPUT  FLUX  IN  THE  STANDARD  OUTPUT 

444  FILE  FORMAT. 

44#  #####  MODIFICATION  #*### 

44#  CRFE  UPDATED  TO  GIVE  MORE  ACCURATE  FLUX  VALUES. 

44#### #####*«#####**#**# **########■#  #»*##**#  44#  #*4##**#»#«# 4###* 

DIMENSION  ENERGY (20), FLX(20),  EPRIME(44> 

DIMENSION  AMASS! 20), RIGID (200), TRANS (200) 

INTEGER  WEATHER 

DATA  ENERQY/30,  ,  40,  ,  SO.  ,  60,  ,  70,  ,  SO,  ,  100.  ,  120.  ,  150.  , 

N  200.  ,  280.',  300,  ,  400.  ,  SOO.  ,  600.  ,  700.  ,  SOO.  ,  1000.  ,  1200.  , 

6  1900.  ,  2000.  ,  2300.  ,  3000.  ,  4000.  ,  5000.  ,  6000.  ,  7000.  ,  SOOO.  / 

DATA  EPR3ME/1.  ,  1.  2,  1.  5,  2.  ,  2,  5,  3.  ,  4.  ,  5.  ,  6.  ,  7.  ,  S.  ,  10,  ,  12,  , 

«<  15.  ,  20  ,  25.  ,  30.  .  40.  ,  30.  ,  60,  ,  70.  ,  SO.  ,  100,  ,  120.  ,  150,  ,  200.  , 

!<  230.  ,  300.  ,  400,  ,  500,  ,  600.  ,  700,  ,  300.  ,  1000.  ,  1200.  ,  1300.  , 

fli  2000.  ,  2900.  ,  3000,  ,  4000.  ,  5000,  ,  6000,  ,  7000.  ,  SOOO.  / 

DATA  AMASS/1.  0079,  4,  0026,  6.  94,  9.  0122,  10.  SI,  12.  Oil, 

*  14.  0067,  13.  9994,  10.  9904,  20.  17,  22,  9890,  24.  305,  26.  9B1S, 

St  28.  0B3 5,  30.  9730,  32.  06,  33,  433,  39.  940,  39.  0983,  40.  OS, 

Si  44.  9339,  47.  9,  30.  9413,  31.  996,  34.  93S,  S3.  047,  30.  9332, 

Si  30.71/ 

OPEN ( UN  IT- 1 , F I LE- ' FLUXES ,  DAT ' ,  STATUS- ' NEW ' ) 

OPEN ( UN  IT-2, F I LE- ' OUTF ILEO.  DAT ' , STATUS- ' NEW  '  > 

444444 444####### #4 #44 #444#44 4 444444444 4444 4444 4444 444 4 4 44# 4# 44 
444  ACCEPT  INPUT  DATA  FROM  COMMAND  FILE: 

44#  YEAR  AND  WEATHER  ARE  DISCUSSED  IN  CRFE. 

444  ORB IT-0  FOR  AIR  PROPAGATION.  -1  IF  GEOMAGNETIC 
444  TRANSMITTANCE  FUNCTION  IS  TO  BE  FOLDED  IN. 

444444# 4*444444444444444 444*4#4444444444444444444444« 444444444 

ACCEPT  4,  YEAR,  WEATHER,  ORBIT 
IF  (ORBIT.  EQ.  1)  THEN 

OPEN ( UNIT-3.  FILE- ' GTRANS.  D AT ' , STATUS- 'OLD  ' , READONLY ) 

DO  J-1,200 

RE AD ( 3 , 20 )  RIOID( J),  TRANS(U) 

END  DO 

10  FORMATUX,  F6.  3,  9X,  FS.  6) 

END  IF 

444444  44444444444444  44444444444444444444  4444444444  44  44444*4*4  4  4 

444  COMPUTE  FLUX  AT  HIGHEST  28  OF  44  ENERGIES.  FOLD  IN 
#4*  GEOMAGNETIC  TRANSMITTANCE  IF  NECEG3ARY.  OUTPUT  A9 
#**  FLUXES.  DAT  (HIGHEST  ENERGY  FIRST). 

C  ##**##***#*#** 4#*###* ##*##*#4#*##**##*444*#*#*#* ##4##*##*##**** 

DO  UG-1,20 
J-29-UQ 
E-ENERGY  <  J ) 

DO  K,»t,  2B 
Z-FLQAT  OS) 

FLX  (K)  -CRFE(  Z,  E,  YEAR,  WEATHER) 

IF  (ORBIT.  EO.  1  )  THEN 
L-l 

RR-AMASS(K) 4S0RT(E4(E+1Q63. ) t/Z/1000. 

DO  WHILE  (RIOID(L) . LT. RR) 

L-L+l 
END  DO 

FLX  (  K)=>TRANS(L )  *FLX  ( K) 
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ENDIF 
END  DO 

FLX  ( 26)  -FI. X  ( 26  )  +FLX  < 27 )  +FLX  ( 23) 

WRITEUi  1)  IFLXIlli  I-t,26> 

1  FORMAT (4<7( 1PEI0. 2)/) ) 

END  DO 

*4*#** ******  #** *44*44*******  ******************  #4**#***#*#****#* 
#4*  COMPUTE  FLUX  AT  ALL  44  ENERGIES.  FOLD  IN  GEOMAGNETIC 
#44  TRANSMITTANCE  IF  NECESSARY.  OUTPUT  AS  OUTFILEO. DAT 
*4#  (LOWEST  ENERGY  FIRST). 

#*  4444**4  *###**##44#4##4  ******  *****  4####**#'*##  **********  ##*###* 

DO  0-1 . 44 

E-EPRIMEI J) 

DO  K-li  20 
Z-FLOAT  <  K ) 

FLX (K )  -CRFE  (  Z,  Ei  YEAR,  WEATHER) 

IF  (ORBIT.  EG,  t>  THEN 
L-l 

RR-AMASB(H ) #BQRT<  E*  (E+ 1863.  i >/Z/1000. 

DO  WHILE  (RIGID(L) . LT. RR  ) 

L-L+l 
END  DO 

FLX (K ) -TRANS (L  >  *FLX (K  > 

ENDIF 
END  DO 

FLX ( 26) -FLX ( 26 ) +FLX <  27  >  +FLX  <  2B  > 

WRITE<2,  1)  (FLX  < I ) ,  1-1. 26) 

END  DO 
END 


FUNCTION  CRFE ( Z,  EZ,  Y,  M) 

##4***4###*4#4#***4*###*##**#*#*#*####*##**##4#44*#4*********** 

4**  THIS  ROUTINE  RETURNS  THE  DIFFERENTIAL  FLUX  TN  PARTICLES/ 

##*  M**2.  STER.  SEC.  MEV/U  IN  THE  INTERPLANETARY  MEDIUM 

*#*  NEAR  EARTH 

#**  FOR  IONS  OF  ATOMIC  NUI4BER-Z 
##*  FOR  ENERGY  (IN  MEV/U)«EZ 
***  FOR  YEAR-Y 

»#*  AND  FOR  WEATHER  CONDITION-M 

44* 

44*  THE  INTERPLANETARY  WEATHER  CONDITIONS  CONSIDERED  ARE: 

*44  M-li  GALACTIC  COSMIC  PAYS  (OCR)  ONLY 

444  M-2i  OCR  +  FULLY-IONIZED  ANOMALOUS  COMPONENT 

*4*  M-3i  OCR  +  90X  WORST  CASE  SOLAR  ACTIVITY 

#**  M«4i  OCR  *  SINOLY-IONI ZED  ANOMALOUS  COMPONENT 

*4*  M-?i  PEAL  ORDINARY  Fl.ARE  FLUX  AND  MEAN  COMPOSITION 

444  M-6,  PEAK  ORDINARY  FLARE  FLUX  AND  WORST-CASE  COMPOSITION 

44*  m«7i  PEAK  WORST-CASE  FLARE  FLUX  AND  MEAN  COMPOSITION 

44*'  M«9i  PEAK  WORST-CASE  FLARE  FLUX  AND  WORST-CASE  COMPOSITION 

*»»  M«9|  PEAK  ANOMALOUS  FLARE  FLUX  AND  MEAN  COMPOSITION 

*4*  M-lOi  PEAK  ANOMALOUS  FLARE  FLUX  AND  WORST-CASE  COMPOSITION 

*44  H-lli  MEAN  ORDINARY  FLARE  FLUX  AND  MEAN  COMPOSITION 

*4*  M«12j  MEAN  ORDINARY  FLARE  FLUX  AND  WORST-CASE  COMPOSITION 

44*  M-10i  MEAN  WQR3T-CASE  FLARE  FLUX  AND  MEAN  COMPOSITION 

*4*  M»14i  MEAN  WORST-CASE  FLARE  FLUX  AND  WORST-CASE  COMPOSITION 

*4*  M«l5i  MEAN  AN0MAL0U3  FLARE  FLUX  AND  MEAN  COMPOSITION 

*44  M-16,  MEAN  ANOMALOUS  FLARE  FLUX  AND  WORST-CASE  COMPOSITION 

*44  Z  IS  the  PARTICLE  ATOMIC  NUMBER 

44*  EZ  IS  THE  ENERGY  IN  MEV/U 
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**«  Y  IS  THE  'r'EAH  OF  THE  . .PECTHA; 

***  1973.  144  «  SOLAR  MIN  ,  1980.  39B  -  SOLAR  MAX 

***  THIS  ROUTINE  IS  BASED  ON  THE  MODEL  IN 
*#*  "COSMIC  RAY  EFFECTS  ON  M ICRCELECTRONICS"  BY 
»*#  J.  H.  ADAMS,  UR,  ,  R.  SILBERBERQ,  AND  C,  H,  TBAQ 
***  NRI.  REPORT  NU.  4506,  AUGUST  25,  1981,  EQUATION  NOS. 

***  AND  TABLE  WOO  REFERRED  TO  BF.LOW  ARE  IN  THE  APPENDIX. 

*******  «*****»**«• **** «* ft*#*#** *#* ******** 4 #4* ***************** 

DIMENSION  AO (3),  ECHO),  B<3> ,  XI  (  6),  X2<3),  CKO),  C2<3> ,  R<92), 

.1,  INDX<92>,  FR(92,  2),  £PX<40>,  FPX  (40> 

***#**«##***■#■** ###**4 *********  * ******* »->«  •»*#*«  i,*#**  **#*»««»*»* 

***  XTSAO  IS  THE  C.  H,  TSAO  CORRECTION  FACTOR  (12/01) 

4* * *********** ************** ************** ************  ********* 

DIMENSION  XTSAO ( 29 ) 

DATA  XTSAO/.  69,  15*.  79,  3*.  OS,  .  47,  3*.  OB.  0*.  47/ 

DATA  IBT/O/i IPTB/0/ 

************4************************************************** 

*##  THESE  ARE  THE  COEFFICIENTS  IN  TABLE  1  OF  THE  REPORT 
************************************  **#******#*#*•»******#******* 

DATA  <  AQ(  I  > ,  I-1.31/-2,  2,-2.  25,  -2.  70/ 

DATA  'EOin,  1-1,31/1  175E5 ,  7,  94E4,  1.  1E5/ 

DATA  (B(I).  1-1.01/2,  73,2.  0,2.  3/ 

DATA  (Xl(I>,  1-1,61/  117,  ,  079,  ,  22,  ,  155,  .  140,  .  117/ 

DATA  (X2(I),  1-1,31/. SO, .BO, .65/ 

DATA  (Cl(I),  1-1,31/6.  52.5,  0,7.0/ 

DATA  (C2(  I  1,  1-1,31/4.  ,  5.  ,  8,  / 

******  ************  I,**#****************#**************#***##**** 

***  ELEMENTAL  RATIOS  FROM  TABLES  2, 3, 4, AND9 
»**#******#*»■*»*+*#**#*■»*#**###*##*****'»«***#****************#* 
DATA  (R(1 1 ,  l-l.  2S1/1.  ,  1.  ,  .  33.  ,  17,  ,  5,  .  023,  1.  .  .  023,  4.  IE-4, 

3.  3E-3,  7.  OE-4,  4.  7E-3.  3.  3E-4,  3,  5E-3,  2,  E-4,  7.  4E-4,  .  07,  .  13, 
it  ,  09,  ,  23. .  05,  .  14,  .  07,  .  14,  .1,1.,  A.  E-3,  4.  0E-2/  .• 

#*#****#*#** »*****#*##**#****#*4#*###*##***#******##*##******** 
***  THE  ELEMENTAL  RATIOS  HAVE  BEEN  EXTENDED  TO  URANIUM  US1NQ 
***  THE  HEAO-3  DATA  INTERPRETED  WITH  CAMERON '3  ABUNDANCES 
*************************************************************** 
DATA  (  R  ( I  ) ,  1-29,921/2.  E-4,  5.  E-4,  4,  E-S,  1.  SE-4,  6.  E-6,  5.  E-5, 

*  a  E-6,  3.  E-5,  3.  E-6.  4.  E-3,  2.  E-6,  2.  E-3,  1,  E-6,  1.  3E-3,  0.  ,  2.  E-6, 

*  4.  4E-7,  Q.  E-6,  1.  E-7,  3.  E-6,  2,  IE-7,  6.  E-6.  3.  4E-7,  3.  E-6,  4.  E-6,  7.  E-6, 

*  4.  E-7,  6.  3E-6,  4.  E-7,  2.  2E-6,  2.  E-7,  2.  E-6,  0.  ,  1.  8E-6,  2.  E-7, 

*  1  SE-6,  2.  E-7,  2.  E-6,  2.  E-7,  1,  2E-6,  2,  E-7,  1.  3E-6,  1.  E-7,  9.  E-7, 

*  9.  E-B,  9.  E-7,  9.  E-9,  1.  E-6,  S.  E-7,  l.  6E-6,  2,  E-7,  1 .  4E-6,  2.  E-7, 

*  8.  E-7,  1 .  3E-7 ,  0.  ,  0.  ,  0.  ,  Q.  .  0,  ,  0.  ,  3.  E-S,  0.  ,  3.  E-3/ 
**************************  ****** ******************************* 
***  INOX  ASSOCIATES  EACH  ELEMENT  IN  COSMIC  RAYS  WITH  A  SET  OF 
***  COEFFICIENTS 

*****************************  *********************************** 
DATA  t  INDX(I),  1*1.  92)/ 1,  15*2,  76*3/ 

**********»***„ ************************************************ 
***  ■  ORDINARY  AND  WORST-CASE  SOLAR  FLARE  ABUNDANCES,  TABLE  6 
****************  ************************************ *********** 
DATA l  (  FR(  I,  J) ,  I-l.  2B1,  J-l,  21/1.  0.  .  022,  0.  0,  0.  0,  0,  0, 

*  1.  6E-4,  3.  SE-S,  3.  2E-4,  0.  0,  3.  IE-5,  1.  6E-a,  4  SE-3,  3.  SE-6, 

*  3.  9E-5,  2.  3E-7.  I,  SE-9,  1.  7E-7.  3.  9E-6,  l.  3E-7,  2.  3E-6,  0.  0,  1.  E-7, 

*  0.  Q,  5.  7E-7,  4.  2E-7,  4.  16.-5,  l.  E-7,  2.  2E-6,  1.  0.  3.  3E-2,  0.  0,  0.  0,  0.  0, 

»  4.  E-4,  l  .  1G-4,  1.  OE-3,  0.  0,  1.  9E-4,  6.  IE-6,  1.  BE-4,  1.  4E-5,  1  6E-4, 

*  1.  IE-6,  a.  4E-3,  a.  E-7,  1  QE-3,  6.  E-7,  1.  E-3,  0.  0,  3.  E-7,  0.  0,  3.  2E-6, 

*  2.  3E-6,  2.  3E-4,  5  5E-7,  1.  2E-3/ 

C  **************************************************************** 

C  ***  THE  ORDINARY  FLARE  ABUNDANCESARE  EXTENDED  TO  URANIUM 
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0  ***  IJSTMO  CAMERON'S  TABLE  OF  GENERAL  ABUNDANCES. 

C  #■#*#*  iHHMHMM*  )l-***4Hh*  ##•###  #+■#**•*****  ###*##*##*##*#****##4*-**#**** 

DATA  (FR<  I,  1  >.  t«29,  92>/2.  E-9,  6.  E-S.  2.  E-9.  3,  E-9,  3.  E-10,  3.  E-9. 

*  4,  e-10,  2.  E-9,  3.  E-10,  1.  E-9,  2.  E-10,  3,  E-10,  4.  E-l  1,  2,  E-10, 

*  0.  ,  9.  E-l  1,  2  E-l  1 ,  6.  E-l  1.  2.  E-l  1 ,  7.  E-l  1, 9  E-l 2,  2.  E-10,  1 .  4E-1 1, 

*  3.  E-10,  6.  E-ll,  2.  7E-10,  2.  E-U,  2.  E-10,  2.  E-l  1.  3.  E-ll,  8.  E-12, 

*  4.  E-l  1,  0.  ,  l.  E- 11,  4,  E-12,  2.  E-l  1,  3.  E-12.  2  E-l  1,  4  E-12,  1.  E-l  l, 

*  2.  E-12,  9.  E-12,  2.  E-l 2,  3.  E-12,  9.  E-13,  1.  E-ll,  2.  E-12,  3  E-ll, 

*  3.  E-ll,  6.  E-ll.  1.  E-ll.  1.  E-ll,  9.  E-12,  1.  E-10,  b.  E-12,  0,  ,  0.  ,  0.  , 

*  0.  ,  0.  .  0.  ,2.  e-12,  0.  ,  1.  2E-12/ 

C  **«  *** ******** ****** ****************  **************** *»***»«  **»* 

C  ***  WORST  CASE  ABUNDANCES  WILL  BE  EXTENDED  TO  URANIUM  BY 

C  ***  SCALING  the  ordinary  abundances  ACCORDING  TO  AN 

C  4*#  ENHANCEMENT  FACTOR  OF  0.222. 

c  *****#*■*#*******##>,  ***************************  ***************** 

DO  300  1-29,  92 
FR  ( 1 ,  2  1-PR  ( 1 ,  1)*0.  22*1 
300  CONTINUE 

C-EZ 

C  *************************************************************** 

C  ***  THE  PARTICLE  MODEL  ONLY  QOEB  DOWN  TO  10  MEV 

C  *************************************************************** 

IF (EZ.  LT.  10.  >  E-10, 

IZ-Z  +  3 

C  ****»»*»**„*  it  it********************  ***************************** 

C  ***  SELECT  THE  MODEL  ELEMENTAL  SPECTRUM 

C  *-.  ,**««*»****  i********** ************************************** 

I'  INDX(IZ) 

r  j-i*2 
'i-ij-i 

c  *************************************************************** 

C  ***  EQS.  4,  9,  AND  6  OF  THE  REPORT 

C  *************************************************************** 

XM-CH  t>*EXP(-Xa<n*(AL0Q10(E>  )**2>-C2!I> 

I  FORMAT  (2X,  El 2.  3,  2X.E12,  3> 

AM-AO( I)*( 1-EXP! -XI ( II J*<AL001O(E) ) »*3 (  I) )  ) 

AX-AO ( t)*< l -EXP ( -X l ( I J) * ( ALOOIO(E) )**B(  I> )  ) 

FMIN-lO**XM 

FMAX-FMIN*(E/EO( 1 ) )»*AX 
FMIN-FMIN*(EXEO( I ) ) **AM 

C  ***********  ,,****  ***********************  ****  »***••***»»**»••*»•* 

C  ***  ADD  FULLY-  [  ONI  ZED  ANOMALOUS  HELIUM  HERE 

C  *************************************************************** 

IPl  CM.  NE.  2) .  OR.  (  IZ,  NE.  2)  )OOTO  10 
IF (E,  LT  200.  1KMIN-.  4 
IF(E.  LT,  300.  ) FMAX-.  03 
10  CONTINUE 

C  *************************************************************** 

C  ***  DO  THE  90V.  WORST-CASE  !  I.  E.  SOLAR  MINIMUM  AND  LOW  ENERGY 

C  ***  STUFF,  EQ13.  11  AND  12) 

C  ************»**#**«*«,»**  *******  **********  *****************  ***** 

1F(M.  NE  3 ) QOTO  16 
IF  (E.  GT.  100.  )  GOTO  13 
IFC I .  GT.  1  )GOTO  1  1 
F-FMIN*(  1  S97*EXP  (  -EX 9.  66  )  +  1 .  A4  > 

GOTO  20 

II  F-FM IN* ( 23.  4*£yP(-E/10.  S4I  +  1 64) 

QOTO  20 

15  F-FM IN  *1 . 64 

QOTO  20 

C  *************************************************************** 
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c  ***  pure  galactic  cosmic  rays  in  year  y.  sqs  i,  a,  and  3 

c  ****************************************************  *********** 

16  AX-.  3* (FMIN-FMAX) 

QX-.  3* IFMIN+FMAX i 
W-.  376*<Y-19S0.  6) 

F«AX*SIN(W)+QX 

C  ******#*****•*************************■**»■**************#*■***+** 

C  ***  DO  LI.  BE,  AND  D.  EC.  7  AND  TABLE  4 

C  **********************************************  *****■!►  *********** 

20  IFK1Z.  LT.  3).  OR.  tIZ.  OT.  3>)Q0T0  30 

IF ( S.  OT.  6000.  ) GOTO  23 
F-.  0142*F 
OOTQ  30 

23  F«F*(.  67#E**<-.  443)) 

C  *************************************************************** 

C  ***  DO  NITROGEN.  EQ,  S 

C  *************************************************************** 

30  IFUZ,  NE.  7) OOTQ  40 

F-F*( .  0064*EXP  <  4*(  ALOO 10  ( E  > -3,  1B>**2>  +  ,  0036* 

5.  EXP  ( 9#(AL0010<E>-.  B>**2>  ) 

C  *************************************************************** 

C  ***  ADD  THE  ANOMALOUS  FEATURE  TO  NITRDQEN,  Ed.  13 

C  ****** ********************************************************* 

IF (M.  NE.  2)  00  TO  45 

FT-(  1.  34E-2)*EXP  ( - < ALOO ( E ) -1 ,  79)**2/.  7> 

IF(FT.  OT.  F)  F-FT 
00  TO  43 

C  *************************************************************** 

C  ***  ADD  THE  ANOMALOUS  FEATURE  TO  NITROOEN,  EQ,  14 

C  *************************************************************** 

40  IF<  (M,  NE,  2).  OR.  <IZ.  NE,  B>  )  OQ  TO  43 

FT- (6  OE-2 ) *EXP  < -  <  ALOQ(E  >  -1,  79)**2A  7)/.  023 
1F(FT.  OT,  F)  F-FT 

C  *************************************************************** 

C  »**  DO  THE  SUB-IRON  NUCLEI,  EOS.  9  AND  10 

C  *************************************************************** 

43  IF  t  I Z.  EQ.  20 )  GOTO  30 

IF  ( (  I Z  .  LT.  17).  QR.  (IZ.  OT.  39 )  >  <30  TO  30 
F-F*<  16.  *<  l,  -EXP  (-.  126*E**.  4 )  )  *E**<-.  33)  > 

30  CRFE-R ( IZ ) *F*XTBAO< IZ ) 

C  *************************************************************** 

C  ***  IF  M.  OT.  4  ADD  THE  CONTRIBUTION  FROM  FLAREB 

C  ***************************************************************** 

IF (M.  OT.  4)  00  TO  100 

C  ******************************************* ******************** 

C  ***  THIS  COMPLETES  THE  COSMIC  RAY  SECTION,  RETURN 

C  *************************************************************** 

RETURN 

C  ************** ft***************************************** ******* 

C  »**  THIS  SECTION  RETURNS  SOLAR  FLARE  ENVIRONMENTS 

C  ***  ' 

C  **#  CONSTRUCT  FLARE  CASE  INDEX,  MU,  AND  HEAVY  IDN  ENRICHMENT 

C  ***  INDEX,  IK 

C  ******  ********  ****#*************************•,*»**#******##***  ** 

100  MK-lM+l )/2-2 

IK-M-2* (MK+l ) 

CRFE-O. O+CRFE 

C  ************************** ************************** *********** 

c  ***  BRANCH  TO  THE  FLARE  CASE 

C  *************************************************************** 

00  TO  (110,  120,  130,  140,  ISO,  160),  MK 
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c  #4+444++**+*  w******##***#****##**##*  #4#+4#*+4+++***+*4+#  ##**•»*+ 

C  **+  ORDINARY  FLARES 

C  *#++*+ #4*+++« 4# *#*4+4*4+##44++ 4+4# #++#### 444*4 *+##+*++*#4+*4 4# 4 

no  if (E.  or,  1000.  i  oa  to  soo 

F-(  1  94E+4  )  *  <  EXP  (  -E/27.  S  1  +  173,  *EXP  <-E/4  )  1 

CRFE-FR <  I Z.  1 1' )  i-F+ORFE 

RETURN 

C  4*444# 44*444 4+44 ####44++**+**++*++#4**+44+#++##+44 4*444**44*44* 

c  4*4  90  PERCENT  WORST-CASE  FLARE9-PEAK  FLUX,  EQ.  22 

C  ##***# 4**#**44#4 44*4 4**#*#*4*44#**##44 444* 44*  »**#*#*4**44# 4***4 

120  1F(E.  OT.  1000  )  00  TO  200 

F« ( 1 .  71E+3)#(EXP  (-E/24.  3) +63.  6#EXP (-E/4) ) 

CRFE»FR(IZ, IK ) fF+CRFE 
RETURN 

C  4#  4444  44*4 4* 4444*4 44 4**4 44*4 4*4444 44444  4 4*44*4444+44 *4* *44 4 444 4 

c  4+4  ANOMALOUSLY  LARGE  FLARES-PEAK  FLUX  EQ.  24 

C  4*44+4 #*+**#4**44**++4++*+**4*****+**+**+#++4+*4+++4+4*4#+*4+*4 

130  PR-SGRT< ( ,  Q01*E1**2+<1.  S6E-314E1 
EX-E+. 01 

PRX-SQRT ( < .  001#EX)**2+( 1.  06E-31*EX> 

DPDE- ( PRX-PR 1 »100 

IF(E.  QT.  130.  )  (SO  TO  131 

F**(9,  3E+9) *DPDE*EXP ( -PR*10,  ) 

QO  TO  132 

131  F-(l.  74E  +  3)4DPDE4PR*4<-9.  ) 

132  CRFE-FRaZ,  IK)#F+CRFE 
RETURN 

C  444444 444*44 444 4 444 *4444444*444*44 #4 4 #44444*44 44 4444444**444444 

C  444  ORDINARY  FLARES-MEAN  FLUX,  EQ,  19 

C  +4*444 444 4 44 *444 4+ *44*4*444444# 4*4 444 4 4444*4*444 +++44+44 4 +44 44+ 

140  IF (E.  OT.  600.  1  00  TO  200 

F«<3.  3E+9)*(EXP  (-E/20.  21+307.  *EXP(-E/3.  1 1 
F-F/3.  2E+3 

CRFE-FRUZ.  I K )  +F+CRFE 
RETURN 

C  4 4+4 4 4 444 444 #44 4 #4 +444+4+4+444 4444 44*4 4 444 4444 4*44 4*4444 44444+ 4 

c  44*  90  PERCENT  WORBT-CASE  FLAREB-MEAN  FLUX,  EQ.  20 

C  444**4 444# 44**44 *44444*4444444444444444***4*4**4*44444*4**4**44 

130  IF(E.  QT.  1000  )  00  TO  200 

F«<7.  6E+9)  *  <  EXP  ( -E/30.  > +  163»EXP  < -E/4.  11 
F-F/3.  2E+S 

CRFE-FR ( IZ, IKl+F+CRFE 
200  RETURN 

C  4444 +4 44444 4#4#44+4++4 44 44 4+4444444444444+44+4444+4444 +4 4+ #4 #4# 

C  #++  ANOMALOUSLY  LARGE  FLARES-MEAN  FLUX,  EQ,  23 

C  4+44 #4 4 +4* 4 4 ##4+44 4*4 4 4 44 4 44 #44 4# 4 44 4+44 +44+ 44# 4 44 +4 44 4 4 44  #  It# 4# 

160  IF(E.OT. 1000. )  QO  TO  200 

F-(  2.  37E+U  >*EXP<  (30.  -El/26.  3) 

F-F/3.  2E+S 

CRFE-FR ( I Z, IK) #F+CRFE 
RETURN 

4*4*4 4 #4*4*4* »*»*+ 44*4444# 4*44 4444 ++**4+* **++##* 4**4 #4 * #44 #44+ # 

#+*  THIS  ENTRY  RETURNS  THE  ORB IT-AVER AQED  TRAPPED  PROTON  FLUX 
#*+  at  THE  SKIN  OF  THE  SPACECRAFT.  THE  FLUX  IS  INTERPOLATED 
4*4  FROM  THE  TABULATION  IN  THE  STASS, DAT  FILE. 

444+444+44 #4##4#44  +  44+  4#4444  #44444+4  #4  4444  44#  *#444*444  4444  44  44+ 

ENTRY  PROTON (EZ) 

IF(  1ST.  EQ.  1)  00  TO  204 
1ST- 1 

OPEN  (UNIT-10,  READONLY.  SHARED,  3TATUS- ' OLD ' ,  FILE-'STASS.  DAT ' ) 

1-1 
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201  READ »  1 0.  202.  END-203)  EPX(1>.FPX<I) 

20?  FORMAT ( lXi  F6.  1,  IXiEIO.  3) 

FPX< I >-795775.  *FPXU> 

I-I  +  l 

oo  ra  201 

203  IPTS-I-l 

204  DO  205  1-2, IPT8 

IF(EZ.  GE.  EP  X  ( I )  )  00  TO  205 
ISAV-t 
00  TO  206 
209  CONTINUE 

PROTON-O. 0 
RETURN 

206  PROTON- <  FPX ( IQAV  > -FP  X  < ISAU-l >) * <  EZ-EPX  < I SAU- 1 M /  <  EPX ( 1SAV ) 
*-EPX< IBAV-1 >  >+FPX ( ISAU-1 ) 

RETURN 

END 
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PROGRAM  RANGER 

*******<■**#•***■***•**********************  ************************* 

#**  THIS  PROGRAM  CONVERTS  A  RANGE.  DAT  FILE  INTO  THE 
***  RANG  IN.  DAT  FILE  NEEDED  IN  MATTER  PROPAGATION. 

»**  THE  RANGE.  DAT  FILE  IS  OBTAINED  USING  TRANSFER. FOR 
***  ON  RANGE-ENERGY  TABLES  GENERATED  BY  J.  ADAMS. 

***  THE  STPOU. OAT  FILE  IS  ALSO  NEEDED  -  IT  IS  GOTTEN 
»**  FROM  STOPPING  P IDLER  TABLES  USING  TRANSFER.  FOR 

*#*-*#*  **#*•»**■***  »*** *********** ***#*44  I************************* 

D I  MENS  ION  C<  6,  44  )  ,  S(4,  44) ,  SS  <  44  > ,  EP  <26,  44)  ,  E(44) 

DIMENSION  ZA<26),  ZB<26>,  i A<26),  IB(26> 

DIMENSION  EG  (  26.  44 ) 

DATA  Z  A/I .  ,  4*2.  ,  £*6.  .  10*8.  .  8*ta.  .  26.  / 

DATA  Z  i3/l .  .2.  ,  4*6.  .  2*8,  ,  10*13.  ,  a*26.  / 

DATA  IA/I.  4*2,  2*3,  10*4,  9*5,  6/ 

DATA  IB/1,  2,  4*3.  2*4,  10*3,  8*6/ 

0PEN(UNIT«3.  FILE- 'RANGE.  DAT ',  READONLY,  STATU9-  'OLD  '  ) 

0PEN(UNIT-2,  FILE-'STPOW.  DAT',  READONLY,  STATUS- 'OLD  ' ) 

OPEN!  UNIT -7,  F ILE-  'R  ANG1N.  DAT',  STATUS- '  NEW  ' ) 

C  **************** *********************************************** 

C  ***  READ  IN  RANGE  DATA  AND  ZERO  OUTPUT  MATRIX 

C  ****** ******************************  *************************** 

DO  10  J»1 ,  44 

READ  < 3,  1  )  £  <  J ) ,  (  C  (  I ,  U  ) ,  1-1,6) 

1  F0RMAKF8.  3,  6(3X,  1PE9.  3>  ) 

DO  I -1.  26 
EP  ( I,  J)-0. 

END  DO 
END  DO 

C  **************************************************** *********** 

C  ***  FIND  ENERGY  CORRESPONDING  CURRENT  PARTICLE  ENERGY  PRIOR 

C  ***  TO  PASSING  THROUGH  1  QM  MATERIAL. 

C  *************************************************************** 

DO  J- 1 , 26 
ZH»ZE(U) 

ZL-ZA(J) 

Z»FLOAT<  J) 

IL-IA(J) 

IR-IB(J) 

DO  K  -1  ■  44 

IF'IH.  EG.  IL)  3<  1,  K>»C(  IL.  K) 

IF  ( IH.  NE  IL)  St  1,  K)-<  (ZH-Z;*C(  IL,  K  >  ♦<  Z-ZL > *C <  IH.  K  >  )  /  (ZH-ZL) 
SS(K)»S(  1.  K> 

FND  DO 

CALL  CU3SPHE,  S,  44) 

DO  K«l,  44 
R  — SS<  K  >  *  1 . 

DO  40  IK-K,  44 

IP(SS(  IK).  OE.  R)  GO  TO  50 
END  DO 

50  IK-IK-l 

EH»E<  IK+1  > 

EL«E<  IK) 

DO  KK-l,  10 

EE— .  5*  <  EL*EH ) 

RP-  K  <  S  <  4,  I K )  *EE+S  ( 0,  1K>)*EE+S<2,  IK))*EE+S(I,  IK) 

I P< RP .  LE.  R  )  EL-EE 
IF(RP,  GT.  R  )  EH-EE 
END  DO 

EP  <J.  K)-.  5* ( EL+EH) 

END  DO 
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END  DO 

****** *4# *44 #4* 4 *4*4 *4 #*+***4 ******** 44**4 ****** 4**4 4**4 44 *4 4 4* 

***  READ  IN  STOPPING  POWER  DATA 

###*##**4#*##*****##*###***#*#***#***#***#**4*#*4*44#***4***##4 

DO  J-l, 44 

REA012,  1)  EiJ),  (C(t.  J  >  I  I-l, 6) 

END  DO 

#4***4  4*44  44  4  *44  4444444#  *44**4  4***4  *4**4 *#**##*##***# 4*4*  4**4  44 

4*4  FIND  RATIO  OF  STOPPING  POWER  PRIOR  TQ  1  CM  MATERIAL,  TO 
*4*  CURRENT  STOPPING  POWER.  THIS  FACTOR  ACCOUNTS  FOR  ENERGY 
*4*  INTERVAL  SPREAD  IN  COMPUTING  DIFFERENTIAL  FLUX. 

*4 4*44 444444 4 44#444444#4*#44**4#*44#4#44#4*4#4 44*444****44 *4*4* 

DO  J»1  , 26 
ZH«ZD(J>**2 
ZL-ZAt J)#*2 
Z-FLQAT< J>**2 
IL-I A( J ) 

IH-IB(J) 

DO  K«  L,  44 

IF  ( 1H.  EG.  IL)  S(l,  K)«CUL,  K> 

IF  (IH,  NE.  IL)  SI1,  K)-<  <ZH-Z)*C(  IL,  k>  +  i  Z-ZL )  #C  ( I H.  k)  >/(  ZH-ZL) 
SS<K)-S( 1,  k> 

END  DO 

CALL  CUBSPLlE,  S,  44) 

DO  K-l,  44 
R-EP1  J,  K) 

DO  IKK-1,44 

T  T  U  u 

IF(tUK).  LE.  R)  GO  TO  ISO 
END  DO 

150  RR«((S(4, IK)4R+3(0, IK))*R+S(2, IK))*R+S(1, IK) 

EG < J,  K)-RRXSS(K) 

END  DO 
END  DO 

************** *#4*#*****4*#**4#*4#**4**44***#***'*#***###*#4**4# 
4*#  OUTPUT  ENERGY  PRIOR  TO  I  GM  OF  AIR 

444444 *444 44 *444*44 44 444 4444 #444*4*4***44*4 4444 4 44 *444 4444 4 4444 

DO  K-l.  44 
KK«*45-K 

WR ITE (7,  2 )  E<K),  (EP(UiK),  U— 1  ■ 26 > 

END  DO 

C  4444*4  44444*4  44  444*444444444#  #444# 44 #44* 444* 4 4 44 444 4444*4 44 4 #44 

C  ***  OUTPUT  RATIO  OF  STOPPING  POWERS 

C  *#*#*# ##«******#**###*#***#*#**##4**#****4*#4#**#**##**#*4**#*4 

DO  K-l,  44 

WR  I T  E  <  7 ,  2 )  E  <  K  ) ,  <EQ(J,  K),  J-l,  26) 

2  FORMAT (4(7<1PE10.2)/)) 

END  00 
END 
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PROGRAM  POWER 

*******#********»****#*#***■*+***********■**************■»*■******* 
***  THIS  PROGRAM  CONVERTS  A  STPOW.  DAT  FILE  INTO  THE 
***  POWIN.  DAT  PILE  NEEDED  IN  HATTER  PROPAGATION.  THE 
***  FILE  SOURCES.  DAT  (AN  ISOTOPE  LIST)  IS  ALSO  REQUIRED. 

***  THE  STPOW.  DAT  FILE  IS  OBTAINED  USING  THE  TRANSFER.  FOR 
***  PROGRAM  ON  STOPPING  POWER  TABLES  CREATED  BY  J.  ADAMS. 
*************************************************************** 
DIMENSION  W<23,  104) , E < 28 > , DE < 2B) ,  C<6,23) 

DIMENSION  12(104). IA(104>, Z1I26),  Z2I26), 11(26), 12(26) 

DATA  Zl/1.  ,  4*2.  ,  2*6.  ,  10*8.  .  B*1B.  ,  26.  / 

DATA  22/1.  .  2.  ,  4*6.  ,  2*8.  ,  10*18.  ,  9*26.  / 

DATA  11/1,4*2,2*3,10*4,3*9,6/ 

DATA  12/1, 2,  4*3,  2*4,  10*5,  0*6/ 

OPEN  ( UN  IT*  1 1  FILE"  'S0URCE3.  DAT ',  READONLY,  STATUS-  'OLD  '  > 
0PEN(UNIT-2,  FILE- 'STPOW.  DAT',  READONLY,  STATUS- 'OLD ' ) 

OPEN  ( UN  IT-4,  FILE-  'POWIN.  DAT',  STATUS- 'NEW  ' ) 

************ *************************************************** 
***  READ  IN  ISOTOPE  LIST 

*************************************************************** 
DO  J-l,  104 

READ  (1,1)  IZ<  J),  IA<  J ) 

FORMAT  (  13,  2X,  13) 

END  DO 

*##***********»********•************«**»•******#•***#**#****»»* 
***  READ  IN  STOPPING  POWER  DATA  IGNORING  16  LOWEST  ENERGIES 
*************************************************************** 
DO  J-l,  16 

READ  (2,2)  01,  QZ,  Q3,  Q4,  QS,  GA,  Q7 
END  DO 
DO  J-l, 2B 

READ  (2,  2)  E<  J),  <Cf  I,  J),  1-1,6) 

FORMAT  ( F0.  1,  1P6E12.  3) 

END  DO 

*************#*******###************************'*************** 

***  DETERMINE  INTERVALS  BETWEEN  ENERGIES 

****** **** ******** ****** ****** ********** ****** ***************** 

DO  JJ-1,27 
J-23-JJ 

DE  <  J  )  "E ( J+ 1 ) -E I J ) 

END  DO 
DE  OS )  -1000. 

****** ************ *********************  *****************  ***■“■*** 
***  INTERPOLATE  STOPPING  POWER  TABLES  TO  ALL  ISOTOPES 
*************************************************************** 
DO  J-l,  104 
I  »IZ ( J) 

IH-Z2 ( I ) **2 
ZL-Z I'D »*2 
Z-FLOAT ( 1 )**2 
IH-I 2 ( I ) 

I L  —  I  1(1) 

DO  100  k-1,  2S 

IF  (IH.EQ.  ID  W(H,  J)-C  (  IL,  K  ) 

IF  (IHNE.  IL)  W(K,  J)»(  (  ZH-Z  )*C(  IL,  K)  +  (  Z-ZL)*C<  IH,  K)  )/<  ZH-ZL) 
W(K,  J)-W(K,  J )  /  I A  (  J ) 

END  DO 

****** «**•****»•****#*•* *************************************** 

***  OUTPUT  STOPPING  POWER  TABLES  TO  POWIN.DAT 

****** ************ ******** **************** ********************* 

DO  JJ- 1 , 20 
J-29-JJ 

WRITE  (4,  4)  E<  J),  DEI  J),  (W(J.  KJ.K-l,  104) 

FORMAT  (FB  1 ,  HX,  FS.  1 .  13 ( /B (  1PE9.  2 )  >  > 

END  DO 
END 
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PROGRAM  APROP 

******  ********  it  ***  **** **************  ********#*#******#****■*•»*** 

***  THIS  PROGRAM  PROPAGATES  A  PARTICLE  FLUX  TAKEN  FROM 
***  FLUXES.  DAT  THROUGH  A  PATHLENOTH  OF  GRAMS 
***  IONIZATION  LOSS  IN  THE  LAST  GRAM  IS  POSTPONED  FOR 
***  THE  PROGRAM  GRAM.  FOR. 

******  ***#*!HHHHHHHt**#***#  *■**#**#****#*♦  #■****•* ***************** 

DIMENSION  F<250, 104) ,  FL1  <240, 104) , FL2( 104), FF<26> 

DIMENSION  FS'104), P< 104),  SPI26),  S( 104) 

COMMON /ALK.JP  IN ( 104),  POUT ( 104) ,  17(  104),  IA( 104) 

COMMQN/BLK/Q ( 104,  104) 

COMMON /CLK/SDEX,  RDEX,  RP,  EMAX,  DX 

******************************  *********‘»************  ***•»*■****■»* 

***  DATA  STATEMENT  DEFINES  TREATMENT  OF  HIGH  ENERGY  PARTICLES 
***  AND  PATHLENOTH  STEP,  DX. 

***  3 DEX “SPECTRAL  INDEX  AT  EMAX 

*•*  RDEX “RANGE  ENERGY  INDEX  AT  EMAX 
***  RP-PROTON  RANGE  IN  AIR  AT  EMAX 

***  EMAX-HIGHEST  ENERGY  TREATED  EXPLICITLY.  BEVOND  THIS  A 
***  POWER  LAW  SPECTRUM  AND  APPROXIMATE  IONIZATION  LOSS 

***  ARE  USED. 

****** ********************************************************** 

DATA  SDEX,  RDEX,  RP,  EMAX,  DX/2.  5,  1.  ,  4  0Q9E3,  9000,  ,  ,  1  / 
***#***************##*******#****#**+****#********•************* 
**#  ACCEPT  PATHLENOTH  IN  GRAMS.  NOTE  DIMENSIONS  OF  F  AND  FL1 
***  MUST  BE  AT  LEAST  F( IORAM, 104) , FL1 < IORAMM, 104) . 
*************************************************************** 
ACCEPT  *,  GRAMS 
IORAM*  INT  (  10.  *GRAMS) 

IORAMM-IGR AM-10 

*************************************************************** 

***  READ  IN  DATA  ON  ISOTOPES  TO  BE  PROPAGATED  FROM  SDURCE3.DAT 
***  AND  STOPPING  POWER  INFORMATION  FROM  P OWN.  DAT. 

****** ******* ***** ********** ****************** ******** ******** » 
OPEN ( UNIT" 1 ,  FILE“ ' S0URCE3.  DAT '<  READONLY,  STATUS-'OLD' ) 

OPEN ( UN IT“4, F ILE“ 'POWIN.  DAT',  READONLY,  STATUS-'OLD' > 

DO  J-l,  104 

READ(  1,  1  )  IZ(J),  I A  ( J ) ,  P  ( J ) 

1  FORMAT!  13,  2X,  13,  32X,  F6,  3) 

DO  K»l,  I  GRAMM 

FL1  (K,  J ) “0. 

END  DO 
END  DO 

**********************************************  *************  »»»* 
***  SET  UP  CROSS  SECTION  MATRIX  (FRAGMENTATION  OF  NUCLEI) 
•*»***.»*•••*»*»»****»***»•***«***«****#•**••**•****•**•*•****#* 
CALL  AIRMAT 

OPEN  ( UNIT-2,  FILE-  'FLUXES,  DAT  ',  READONLY,  STATUS-'OLD'  ) 

DO  1-1. 29 

****** ******************** ************************************* 

***■  READ  IN  ELEMENTAL  ABUNDANCES  FROM  FLUXES, DAT  AND  CONVERT 
***  INTO  I3QT0PIC  ABUNDANCES  USING  APPROXIMATE  ISOTOPIC 
***  FRACTIONS  P ( 1 04 )  FROM  S0URCE3.DAT. 
**************************************************************  » 
READ  (  2,  2  )  (SP(K),  K-l,  26) 

2  FORMAT <4(7( 1PE10.  2)/) ) 

DO  K-l,  104 

IZZ^IZIK! 

S I  K )  »£3P  ( I ZZ )  *P  < K ) 

FL2  (K )  -S  ( K ) 

DO  J-l . ICHAM 
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F<J.  K>«0. 

END  00 
END  DO 

4##44#4*###4**#4###########4#4*4#4###4##44*###*##*444#4*4*###** 

**#  SET  UP  IONIZATION  LOSS  VECTORS 

44*444 44**## 444444444*44*44*44 4444 4* 4**##44# 44*44**4 444 4 4# 4 *4*4 

CALL  ILQS8 

4444 44 444**4 44 4* 44 #4 4* 4*4# 4*44 *44# 44 44 44444* 4* 4# *4*4 4 444 *4 44 44* 

*4*  FRAGMENTATION  AND  IONIZATION  LOSS  FOR  ALL  BUT  LAST  GRAM 

#44 4 44 4*4444 44*4 444**4 4 4 *4 #44* 4#*###4 *4*44 44 44 *444 4* #44#  4*4##*# 

DO  J*l.  IORAMM 
DO  KA«1. 104 

F8<KA)»P!N<KA)*FL1 ( J. KA ) +POUT (KA > *FL2 < KA > 

END  DO 
DO  KA-1. 104 
DO  KD-1, 104 

F<  J.  KA)-F(  J,  KA )  +G<  HAi  KB )  *FS  ( KB  ) 

END  DO 

FL2<KA)«F<  J,  KA) 

END  DO 
END  DO 

#####* 4*##*#*#*#4#**#####*#4**###*44*4###*4**4####**###*#*##4#4 

44#  SET  UP  BOUNDARY  CONDITIONS  FOR  NEXT  LOWER  ENERGY 

444444 4 444444444 44444444 4444444444444444444 4444#444###4#4 4 4444# 

DO  J-l,  IQRAhM 
DO  KA" 1  j  104 

FL1IJ,  KA)«F<J.  KA) 

END  DO 
END  DO 

#44 444 #444444444 44 4444444*444*44 444 4 444444 4444444444*4444444444 

**4  PERFORM  FRAGMENTATION  ONLY  FOR  LAST  ORAM 

4444 4 4 4444 444 44444444444 44444*444444 *4 44444*44 4 444 #444 4*444444  4 

DO  J-IORAMM+l,  IQRAM 
DO  KA-l.  104 
DO  Kfl-1.  104 

F<  Ji  KA)»F<  Ji  KA)+(HKA,  KB)#F<J-l,  KB  ) 

END  DO 
END  DO 
END  DO 

444*4# 4444 444* 4444 4444 4444444 4 444 444 #4444 444 4*4444 44 4*4 44 4  44444 

*4#  PUT  SUM  OVER  ELEMENTS  INTO  FF(26> 

444444 #44444444 44444444444444444 444444444444444 44 444444444 4444 4 

DO  JZ-1.  26 
FF  <  JZ ) “0. 

END  DO 
DO  K-l,  104 
IZZ-IZ<K> 

FF ( IZZ  >-FF< IZZ >+F< I  ORAM. K> 

END  DQ 

44*44# 4*44*44444444444 4#44#*4444444*44444#44444444*4#4 4444 44*4 4 

444-  OUTPUT  RESULTS  FOR  CURRENT  ENERGY 

44*444 4444*44*44 444*4# 44444*4***444444444444 44444444 4*4444 4444 4 

WRITE(IO.S)  (FF(R), K-l. 26) 

FORMAT)  4  (  7  (  1PEIQ.  2)/)  ) 

#*#*#* 4*#### 4*4# 4 *#***#44*4** 4 4*****4##» 4*44 4 #**#***#**#*#*##* 4 

*#*  PROCEED  TO  NEXT  OF  2B  ENEROIES 

##**##  44*4*4 4**4 4**4**##*##**# 44# 4 4**4**#* #4* 4 4 #4# 4# 44*4 ******* 

END  DO 

END 
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SUBROUTINE  AIRMAT 

4*4*44#  4*4 ff**#**#*#**##*#*####*#**#**###***#** ********  *444-4* **•* 
#**  THIS  SUDROUTINE  SETS  UP  THE  CROSS  SECTION  MATRICES 

#*#444#  ft###**#**####  44  44  4  444444444444444444444  4  4444444*444  4444  4 

COMMON/ ALK/PIN( 104),  PdUT( 104 > .  1 2 ( 104 >  , I A ( 104 ) 

CQMMQN/ELK/Q<  104,  104) 

COMMON /CLK/SDEX,  RDEX,  RP,  EMAX,  DX 
DIMENSION  10(104) 

44#444#444*4444#44444444#*44444444444444444444*444#4*4#44444#.*4 

*4*  READ  IN  PARTIAL  CROSS  SECTION  DATA  PROM  0.  DAT 

4***4#  *444*4 44 4444444*444444*4*4 #444444*444*44*444*4*4 444* *44 44 

0PEN(UNIT»3,  FILE*  '0.  DAT',  READONLY,  STATUS- 'OLD  '  ) 

DO  J-l ,  104 

READ  <  3,  1  )  (  I Q  (  I  ),  [  =  1,  104) 

1  F0RMAT(1X/,  6(2014/) > 

DD  J-l!  104 

Q!  I ,  J ) -FLOAT < I Q < I )  ) 

END  DO 
END  DO 

****** ******44*4***# ********************************4*##***#*** 

#**  REPLACE  CROSS  SECTIONS  WITH  RECIPROCAL  MEAN  FREE  PATHS 

*##**4  **************************** 4***#*****##*##**#4**#******4 

FA-5.  97SE-4/14.  42 
DO  J-l,  104 
DO  I-i,  104 

0<  I,  J)  -FA*DX#G(  I,  J) 

END  DO 
END  DO 

##**#***#***#»****************#**#***#»*************#*******.*** 
***  INCORPORATE  PROTON-NUCLEUS  TOTAL  INELASTIC  CROSS  SECTION 
*#  **************************************************************** 
C1-.  79*44.  9*14.  **.  7 
C2-.  21  *44.  9*16.  *#,  7 
Q<  1  <  1)-1.  *0(1,  1 )  -FA#DX*( C 1+C2) 

*###*#*#**#*****##**#**••#•******#*****#»*#*****#»#«****«***»*# 

***  INCORPORATE  NUCLEUS-NUCLEUS  TOTAL  INELASTIC  CROSS  SECTION 

****************#*##*#***#*##***#**#***#*********#****4*#**#*** 

DO  J-2,  104 

C1-.  79* (  IA<  J)**.  333+14.  **.  333-.  41**2 
C2-.  21*  <  I A  (  J)**.  333+16.  #*.  333-.  4)**2 
Q(  J,  J)-l.  +Q  <  J,  J )  -49  SS#FA*DX  *(  C 1 +C2  ) 

END  DO 
RETURN 
END 


SUBROUTINE  ILC3S 

**#*#**4**4*  44*#*** #*###****4#*** #4*4**##***#**# ***#44*4** »«*## 

44*-  THIS  SUBROUTINE  SETS  UP  THE  IONIZATION  LOSS  VECTORS 

44*44* #*##*4**4 4*4 ##**#4** 4 *4# *4*4 4*4*4 4 4**4***44*4444 44*4 4**44 

COMMON /ALK/P  INI  104),  POUT(  104),  IZd04),  IA(104) 

COMMON  /CLK/SDEX ,  RDEX,  RP,  EMAX,  DX 
DIMENSION  WL(  l')4),W(  104) 

C  *4*44*  #44*44 *444 ***4444444*4 4 4 *4 4 4 4*4 44*44**4* *44*44 4*44*4 4**4 4 

C  ***  READ  IN  STOPPING  POWER  DATA  FROM  POWIN.  DAT  (FOR  CURRENT 

C  +*#  ENERGY  ONLY). 

C  *#*44*  44**44 4 4*4 4444*444 444444 4444 444 44444*444 4**444  44*44*44**4 

READ  <  4,1)  E,  OE,  (Wl  J),  J-l,  104) 
t  FORMAT  (FS.  3,  5X,  FB.  3,  13</9<  1PE9.  2)  >  ) 
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************ *********#****+***+*** ****** *********************** 
***  CALCULATE  IONIZATION  LOSS  VECTORS  USINO  HIQH-ENERQY 
*##  APPROXIMATION  IF  ENERGY  ■  EMAX,  QR  DERIVATIVE 
*»*  APPROXIMATION  AT  OTHER  ENERGIES. 
*************************************************************** 
DEX-USDEX-i.  > *RDEX  +  i .  )/RP*DX 
DER-1.  /DE 
EM-EMAX-.  OS 
IF  <E,  LE,  EM)  THEN 
DO  J»i.  104 

PINCJ>«UL< J)*DER*DX 
POUT ( J ) *1 .  -W< J)*DER*DX 
END  DO 
ELSE 

DO  J-l,  104 
PIN(J>«0. 

POUT! U )“1 .  -FLOAT < IZ( J) ) **2*DEX/ FLOAT < IA(U> > 

END  DO 
ENDIF 

*************************************************************** 

***  STORE  STOPPING  POWERS  FOR  CURRENT  ENERGY  IN  ML (104)  FOR 
***  USE  AS  'PREVIOUS'  ST0PP1N0  POWER  NEXT  TIME  AROUND. 
************ *************************************************** 
DO  d«l.  104 
WL(U)-W<  J) 

END  DO 
RETURN 
END 


J 
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PROGRAM  ORAM 

*4 4444 444444 444 4444444 *444 44444444 #44 4444 44444 444 444444444 *44 

#**  THIS  PROGRAM  CONVERTS  PROPAGATES  THE  DIFFERENTIAL 
##*  ENERGY  FLUX  THROUGH  1  QM/CM442  MATERIAL  TAKING  INTO 
444  ACCOUNT  ONLY  IONIZATION  L0S3.  IT  ALSO  EXTENDS  THE 
***  MINIMUM  ENERGY  DOWN  TO  1  MEW AMU. 

4444444444«4444444*44*#4*444444444*4#44**4444444#444»444*44#4 

DIMENSION  FINl26> 20),  FQUT<26,  44) ,  Et 44 ) , EP (26, 44) 

DIMENSION  F  <  4,  2B ) , EN< 20 ) ,  EO<26, 44) 

444444  44444444444444444444.4444« 4444444444444444444444*4444*4444 

444  READ  IN  FLUX  FROM  FOROIO. DAT 

t,44444444**4#*44444444*44444#44*4444*444444444444444444444**444 

DO  1 1* 1 <  2S 
I-29-II 

READ!  10,  1  )  <F1N(  J,  I  ) ,  J-i,  2t>) 

1  FORMAT <  4  (  7  <  IPE10.  2.)/)  ) 

END  DO 

4#4#444444#*44*444*#**##4####**44*4#********4444#4###4*444*4444 

**«  READ  IN  TABLES  OF  ENERGY  PRIOR  TO  1  GM  AND  ENERGY 
444  INTERVAL  SPREAD  FROM  RANGIN.DAT. 

#44##444444*44*444#4##*#4##4#44*4#444'**4*4*44#4#**44#4*44#*4444 

0PSN<UNIT-7,  FILE"  'RAN01N.  DAT READONLY,  BTATUS- 'OLD ' ) 

DO  I"l,  44 

READ (7,  1 )  E(I),  <EP<J.  I),  J-I.2L) 

END  DO 
DO  1-1,  44 

READ <7.  1 )  Ed),  <EQ<J,  I),  J«1,2S> 

END  DO 
DO  1-1, 2S 

ENd)-Ed  +  16> 

END  DO 
DO  J-I,  2b 
DO  K-I,  4 
F<K,  J)-0. 

END  DO 
DO  I-l,  29 

Fd,  I )  «FIN(  J,  I) 

END  DO 

CALL  CUB3PUEN.  F,  20) 

DO  1*1,4 

F(  I,  20 )-F(  I,  27) 

END  DO 
DO  I  — 1 , 44 
ER-EP  ( J,  I) 

DO  KK-1, 20 
K-29-KK 

IF(EN<K>.  LT.  ER)  GO  TO  30 
END  DO 

FOUT(  J,  I  >  =» <  <  <F<4,  K)4ER+F<3,  K>  >*ER+F<2,  K)  >*ER*F<  1,  K)  )* 

!<  EQ(J,  I) 

END  DO 
END  DO 

444*44 444444 4 4*444444444444 444 44*444 4 #4*44 *44*4444 44 4444 *4 4*4 

444  OUTPUT  DIFFERENTIAL  FLUX  IN  FINAL  FORM 

444444  444444 444444444*44*4 44**4444*4 444*444444 44  44  4444 44 44 444 

DO  J«I,  44 

WRITE120,  1)  (FOUTU.  J),  I-l,  26) 

END  DO 
END 
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PROGRAM  ATMOS 

*************************************************************** 
***  THIS  PROQRAM  CONVERTS  DIRECTIONAL  FLUXES  < "QUT2" )  AT 
***  DEPTHS  OF  0,  5,  13.  30,  30,  73-  AND  100  QM/CM**2  IN 

***  THE  ATMOSPHERE  AND  INTEORATES  THEM  OVER  ALL  ANGLES  TO 
***  SET  THE  TOTAL  FLUX.  DOTH  GEOMAGNETIC  CUTOFF  AND  THE 
***  INCREASED  GRAMMAGE  AT  NON-VERTICAL  DIRECTIONS  ARE 
*#-»  INCLUDED. 

***  INPUT  FILES  ARE  NAMED  AFILE1  THROUGH  AFILE7, 

******  ************************  ************  ********************* 

DIMENSION  F<26, 44) 

CQNHON/ALK/E  (  44 ) ,  C(7,  26,  441 

*****  ********************************************************** 
***  ENTER  ALTITUDE  AND  CUTOFF,  CALCULATE  GRAMMAGE. 
************************ *************************************** 
ACCEPT  *,  ALTITUDE,  CUTOFF 
CUTOFF  "CUTOFF  / 1 0, 

CALL  ALTGRM< ALTITUDE,  GRAMS! 

*************************************************************** 
***  ACCEPT  FLUX  DATA 

****************************** ********************************* 

0PEN<UNIT«21,  FILE-'AFILEl.  DAT',  STATUS* 'OLD ' , READONLY) 
0PEN1UNIT-22,  FILE- ' AFILE2,  DAT'.  STATUS" 'OLD  ', READONLY) 

OP EN< UNIT-23,  FILE- 'AFILE3.  DAT',  STATUS- 'OLD  ' ,  READONLY) 
0PEN(UNIT«24, FILE- ' AFILE4.  DAT',  STATUS- 'OLD ', READONLY) 
0PEN<UNIT-23,  FILE- '  AF ILE3,  DAT ',  STATUS* 'OLD ',  READONLY > 
0PEN(UNIT-2&,  FILE-'AFIL£6.  DAT ',  STATUS- 'OLD ', READONLY ) 
0PEN!UNIT«27,  FILE- ' AFILE7.  DAT ',  STATUS- 'OLD ' , READONLY) 

DO  L-l, 7 
LL-L+20 
DO  K- 1,  44 

READ!  LL,  402 )  (C  ( L,  I,  IX  1 ,  1-1 .  26) 

END  DO 
END  DO 

CLOSE! UNIT-21 ) 

CLOSE! UNIT-22) 

CLOSE! UNIT-23) 

CLOSE ( UNIT-24 ) 

CLOSE! UNIT-23) 

CLOSE! UNIT-26) 

CLOSE! UNIT-27) 

C  ************************************************************** 

C  ***  MAIN  ROUTINE 

C  **** ********************************************************** 

DO  1-1,9 

THETA-FLOAT! I- 1)#.  174333 
GM-0RAM8/C0B! THETA) 

DO  J-l.  12 

IF  ( (  I.  EG.  1.  AND.  J.  EG.  1  ).  OR.  (  I,  OT.  I  )  )  THEN 
PHI-FLOAT! J-l)*. 323399 
CUT-CUTFACT(CUTOFF, THETA, PHI ) 

IF  (I.  OT.  I  )  THEN 

AREA-  S209S8*! COS (THETA-.  087267 1 -COS ( THETA*.  097267) ) 
ELSE 

AREA-. 023909 
END  IF 
DO  K-l, 26 

RESIDUAL-RANGE!  1000.  *CUTOFF,  IXJ-GM 
DO  L-44,  1,-1 

IF ( RANGE! E (L> , K!  OT.  RESIDUAL)  THEN 
P  (K,  L)»P  (K,  D+FLUX  (GM,  IX,  L  )  *AREA 
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ENDIF 
END  DO 
END  DO 
ENDIF 
END  DO 
END  DO 
DO  L-l.  44 

WR I TE  (60,  402 )  <  P  ( K.  L ) .  K« 1 ,  26 ) 
END  DO 

402  FORMAT  (4<7<  IPE10.  2)/ )  ) 

END 


SUBROUTINE  ALTQRMl  X.  0 ) 

C  ***#***#*****■***•»■#■#•**•»•#■*»■*■)>*•»•*■**•»•»*#■**#•»■**#*■»*#•*##*-»*■»***.*■»+*#* 

C  ***  THIS  SUBROUTINE  CALCULATES  THE  QRAMMAQE  CORRESPONDING  TO 

C  #**  A  GIVEN  VERTICAL  ALTITUDE  IN  KILDFEET. 

C  #**#■»**■»** #*■**#*■»#***■»**#*#*■»*•♦*•»**#****#*#'**■»*•*##*■»*#  *■»****##* 

ZOM-.  04S34-1.  17E-9*AG9<  X-10S  >*#3.53 
0-1033.  *EXP(-ZQM*X) 

RETURN 

END 


FUNCTION  RANOE(ENER0V.  IZ) 

*#*«**  4##***'il-*»***'IHMM#***##**-»*-»*#***##**#***«'***.4**  *♦*##**#** 

***  THIS  FUNCTION  CALCULATES  THE  RANGE  OF  CHARGED  PARTICLES 
***  IN  AIR  (79V.  N/21V.  0)  USING  THE  BETHE-3L0CH  FORMULA.  THE 
***  CHARGE/MA93  RATIO  13  TAKEN  AS  2.  ALL  ENERGIES  IN  MEV/NUC. 

**.#**»«**<***«*.i,j**«***.*#«HH(**i****m**»m****tM***4  4M» 

DIMENSION  XUO),  A<10).  AMASS<26! 

DATA  X/-.  97390633,  -.  B6S06337,  67940937,  43339339, 

!<  1 4907434,  .  1 4007434,  .  43339339,  .  67940957, 

h  .  06506337,  .  97390653/ 

DATA  A/.  06667 1 34,  .  14943133,  .  21900636,  .  26926672, 

?<  29  3  32  422,  .  29352422,  .  26926672,  .  21900636, 

V  .  14943135.  .  06667134/ 

DATA  AMASS/ 1.  0079,  4.  0026,  6.  94,  9.  01210,  10.  31,  12.  Oil, 

4»  14.  0067,  13.  9994,  10.  9904,  20.  17,  22.  99977,  24.  305,  26.  90134, 

>A  20.  0055.  30.  974,  32.  06,  35,  453,  39.  940,  39.  09B3,  40.  00,  44.  9559, 

!,  47.  90,  50  9413,  51.  996,  54.  930.  33.  847/ 

»*#***  ***•*■*«*  **»*«.»#*«•*##•*****»•*»***«*#  #  **»»**»  »**»#***«*»# 

###  CALCULATE  RANGE 

#**#**  •*■**#*****■*#  4*#* 

S0-10. 

R0-.  145 
PARITY-1. 

IF(  ENERGY.  LT.  EO)  THEN 
Y-EO 

EO-ENERGY 
ENEROY-Y 
PARITY— 1. 

ENDIF 

SLOPE- .  3*( ENERGY-EO) 

YINTERCEPT-.  5 M ENERGY  +E0 ) 

RANOE-RO 

IF  (ENERGY.  NE.  EO)  THEN 
DO  J-l,  10 

E-SLQPE#X ( J ) +Y INTERCEPT 
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RANGE-RANOE+AI J) »PAR ITY/3TP0W< E ) *SL0PE 
END  DO 
ENDIF 

RANGE-AMAS31  I Z )  »R ANGE/FLOAT  <  1 1 )  **2 

RETURN 

END 


FUNCTION  STPOW(E) 

C  *##•*■**  *#****•!***#■*»•«■*#**#•»*##***#»**#■****#***•**■«******■*#*■»•»#■)*•*** 

C  **#  THIS  FUNCTION  CALCULATES  THE  STOPPING  POUER  OF  PROTONS  IN 

C  ***  AIR  US.  ENERQY  IN  MEV. 

C  **##■*■#  ##*#**•**«■*+#***#■«*■*'**#  *•*■*#*■*****■****■»*■»#■»■****#****** 

GAMMA-1.  +E/931.  5 
0-1 .  -1. /0AMMA*#a 

X1-,  1S344*(  ALC1Q  (  12460,  *D/<  1.  -0  )  >-0>/B 

X2-.  1 3334*1  ALOO (  10376,  *B/(t.  -OM-fll/B 

BTPQW-. 7670S»X1+. 23292#X2 

RETURN 

END 


FUNCTION  CUTFACT (CUTOFF,  THETA,  PHI > 
****************#*************#+***********#**********#*#****** 
***  THIS  FUNCTION  CALCULATES  THE  OEOMAGNETIC  CUTOFF  AT  ZENITH 
*#*  ANGLE,  THETA,  AND  AZIMUTHAL  ANGLE,  PHI.  A  DIPOLE  FIELD 
#**  IS  ASSUMED,  CALCULATION  AT  20  KM, 
***************************************************** ********** 
FACT-39.  6/(1,  +20,  /6371.  )  #*2 
COQL-14.  PCUTOFF /  FACT )  **,  S3 
C0SG-C0S1 PHI )*SIN( THETA) 

A«(  1.  +9GRT1  1.  •■COSG#C03L**3)  )  **2 

CUTF ACT-4.  #CUTOFF/A 

RETURN 

END 


FUNCTION  FLUX (GM. IZ, IE) 

**#****#*•**■**•**##*#**•***•»*****##.*■*#*#■**»#**  i******************* 

***  THIS  FUNCTION  INTERPOLATES  THE  FUNCTION  C(7,26,  44) 

**#  TO  QIVE  ITS  VALUE  AT  ORAMMAQE,  OM. 
**#■*■**»•**»•»****■**■••**#■»***********■*■»***##***■*•*■**■♦■***■*•»■*■*■****■»**•* 
DIMENSION  0(7) 

COMMON / ALK / E ( 44 ) , C ( 7 , 26, 44) 

DATA  E/1.  ,  l,  2,  1,  3,  2.  ,  2,  3,  3.  .  4.  .  3.  ,  6.  ,  7  ,  B.  ,  10  .  12.  ,  13.  , 

?<  20.  ,  23.  ,  30.  ,  40.  ,  30.  ,  60.  .  70,  .  SO,  ,  100,  ,  120.  ,  130.  ,  200.  ,  230.  , 

fj  300.  ,  400.  ,  300.  ,  600.  ,  700.  ,  BOO.  ,  1000.  ,  1200.  ,  1300.  ,  2000.  , 

5  230Q.  ,  3000.  ,  4000.  ,  SOOO.  ,  6000.  ,  7000.  ,  8000.  / 

DATA  D/0.  ,  3.  ,  15  .  30.  ,  SO.  ,  73.  .  100,  / 

DO  J-l,  6 

IF  (GM.  GE.  D(U>)  ,JX«U 
END  OQ 

IF  (C(UX.  IZ,  IE).  LE.  0.  .  OR.  C(UX  +  t,  IZ,  IE).  LE.  0.  >  THEN 
FLUX-0. 

ELBE 

XX-C(UX,  IZ,  IE)  /C  ( JX  +  t  i  IZ,  IE) 

X«AL00(XX)/(D(JX+1)-D(JX) ) 

IF  (X*(OM-D(JX)  ).  LT.  13.  )  THEN 
FLUX-CUX.  IZ,  IE)  -EXP  ( -X*  (GM-Dl  UX  )  )  > 

ELSE 

FLUX-0. 

ENDIF 

ENDIF 

RETURN 

END 
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PROGRAM  LETINTT 

***************************  ********************  *4**4*  *********** 

***  THIS  PROGRAM  CALCULATE3  THE  INTEGRAL  LET  SPECTRUM  SUMMED 
***  OVER  ALL  ELEMENTS  FROM  A  FLUX  (OUTF1LE,  DAT)  FILE. 

***  INPUT  FILE  NAMED  FQRO 10.  FOR 

***  OUTPUT  FILE  NAMED  FOR02O. FOR  (USE  PLOTT  TQ  PLOT) 
************** ********** *************************************** 

DIMENSION  C<6,44),  ZZ(2S).  tZ(£B),S(4,  44),  XE(500>,  DXE(SOO) 
DIMENSION  E(44),  H<4.  44),  G(26,  44),  A (27,  4,  500) 

DATA  Et,  E2/1.  ,  3000,  / 

DATA  22/1.  ,  2*2.  ,  0*6.  ,  7*9,  ,  3*18.  ,  7*26.  / 

DATA  12/1,  2*2,  3*3,  7*4,  3*5,  7*6/ 

DO  J-l . 4 

*************************************************************** 
***  INITIALIZE  SPLINE  PARAMETER  MATRICES 
********************************* ****#***##****##******* ******* 
DO  K« l, *44 
S <  J,  K )  -0, 

H(  J,  K>-0, 

END  DO 

*************************************************************** 

***  INITIALIZE  OUTPUT  LET  MATRIX 

***  K  -  CHARGE  (FOR  K-27,  SUM  OVER  ALL  CHARGES) 

»**  J  -  SPECTRUM  TVPE  (1-INTEGRAl.  LET) 

***  L  -  DATA  POINT 

****** *************************************************** ****** 
DO  K-1,27 
DO  L-l,  500 
A (K,  U,  L)-0. 

END  DO 
END  DO 
END  DO 

**************************************************************** 
***  READ  SILICON  STOPPING  POWER 

******  *******.»4*********  ********************  fl.****  #«•■***■*****■**** 

OPEN  (UNIT-8,  FILE-'BTPOW.  DAT ',  READONLY,  STATUS- 'OLD  '  ) 

DO  J-l,  44 

READ  <  0,  t )  E(J),  (C<  I,  J) ,  I  - 1 ,  6) 

FORMAT ( F8.  3,  60X,  1PE9.  3)  ) 

END  DO 

***************************************«.,  ,,„**#***»#******#*#*»* 

***  READ  IN  FLUX  DATA  FROM  F0R020. DAT 

********************  **************  **4*********  ******  »***«•***«** 

DO  I -1,44 

READ  (20,  2)  (0<J,  I),  J-l,  26) 

F0RMAT(4(7(  1PS10.  2)/  ) ) 

END  DO 

*********************************************************  ******* 
***  DETERMINE  ENERGY  AND  LET  BIN  WIDTH  PARAMETERS 
»**•*• ************************************************** ******* 
Fl-( E2/E1 ) **. 002 
F£-(  1.  -1.  /FI  )  *E1 
FS-  5*  (  1.  +1.  /FI  )  *E  1 
01-10.  **.  01 
02-1.  -1.  /Ot 
DO  1-1 ,  500 

XE ( I ) *F3»F 1  *  *  I 
DXE(  I  )»F2*FU*I 
END  DO 

****** 4 ********* *** *************************** 4**444**** ******* 

***  PERFORM  CALCULATIONS  ON  EACH  ELEMENT  FIRST 
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C  *************************************  *********#****#■********  +  #** 

DO  J-l,  26 

C  #**#»#******  I******************  **************  ***###*#****#*■*#*** 

C  ***  ESTABLISH  FLUX  POINTS  TO  BE  SPLINE  INTERPOLATED 

C  ***********************************************  ***************** 

BO  1-1,44 

H(  1,  I  >-3<  J,  I  ) 

END  DO 

CALL  CUBSPL  ( E,  H,  44) 

C  *******#*****************************************»******#*■„■**** 

o  ***  ESTABLISH  STOPPING  POWER  POINTS  TO  DE  SPLINE  INTERPOLATED 

c  **********************************  ****************************** 

Z 1 -FLOAT ( J ) 

Z2-ZZ ( J ) 

I2Z«I2(J) 

DO  1-1,44 

AA-Wd,  +.  001 073*>27*E  ( I )  >  **S 
13-SQRT  <  1 ,  -AA) 

Z1E-Z1*<  1.  -EXP<-130,  *0/21**.  6667)  ) 

Z2E-Z2* ( 1.  -EXP  <  “130,  *B/Z2**.  6667)  ) 

3(1,  I )  - (  Z  IE  / Z2E)#*2#C  ( I ZZ,  I) 

END  DO 

CALL  CUBSPL(E,  3,  44) 

C  **************************************************************** 

C  ***  FOR  EACH  BIN  USE  FLUX  AND  STOPPINO  POWER  INTERPOLATIONS 

C  ***  (PARAMETERS  CONTAINED  IN  H  AND  S,  RESPECTIVELY  >  TO  QET 

C  ***  INTEGRAL  AND  DIFFERENTIAL  LET  AND  ENEROY  SPECTRA. 

C  *************************************************************** 

1DEX-1 
MDEX-1 
DO  II-L,  300 
I-301-I I 
X-XE(I) 

DX-DXE  < I ) 

DO  KK-IDEX,  44 

IF(E<K) .  LE.  X)  00  TO  40 
END  DO 

40  F- ( (H ( 4,  K )  *X+H  ( 3,  K ) )*X+H(2,  K) > *X+H( 1 , K> 

1DEX-PK 
DO  kk-MDEX,44 
K-49-KP 

IF(E(k).  LE.  X  )  GO  TO  60 
END  DO 

60  SP-< (S(4,  K)*X+8<0.  K) )*X+S(2, K) >*X*B( 1, K) 

MDEX-Kk 

A(  J.  4,  I  )«A(  J,  4,  I  )+F 
DO  k-1,  I 

A  <U,  3,  K  >  —A <  J ,  3,  K  )  +F*DX 
END  DO 

•  L-  l  +  INT  <  100.  *ALQO  10  (  SP  )  ) 

DL»G2*01**L 

0-F*DX/DL 

<iP-F*DX 

A(  J,  2,  L)-A(J,  2,  L)*Q 
Ai27,  2,  U»A<  27,  2,  L)+Q 
DO  K-l,  L 

A<  J,  l,  K>»A(  J,  1,  Kl+QP 
A  (27,  1,  K)=*A(B7,  1,  P)+QP 
END  DO 
END  DO 
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****** 4 *44**44***  **************  **************  *************  *444* 

***  ADD  HIGHEST  ENERGY  COSMIC  RAYS  TO  LET  SPECTRA 

*****#*##*4#444#4**4*##*4*4  4*4***4  *****4#*##**  ******4444*  **4**4 

EP=*<2t/22)**C»C (122,44) 

L«l* INTI  100.  #AL0G10(SP ) ! 

DL“024G1*#L 
G-G(  J.  44>*E2/1.  5/DL 
QP-OtU.  44)*E2/1.  S 
A(  J,  Si  L)-A(Ji  2.  L )  +G 
A  (  27  •  2.  L)«A<27,  2.  L)+Q 
DO  K«l.  L 

A<  J.  t,  K)»A<  J,  1,  Kl+QP 
A <  27,  1,  K)=>A<27,  1,  K)+GP 
END  DO 
DO  K»t,  SOO 

A  (  J,  3i  K )  “A  ( Ji  3i  K  )  +QP 
END  DO  • 

END  DO 

C  #***44***#*#*#*#*#44***#*******4*#*#*##+44**44»**4**4#***4  4***4 

C  ***  DETERMINE  LET  VALUES  AND  INTEGRAL  LET  SPECTRUM  VALUES 

C  ***  AND  OUTPUT  TO  FQH030.  DAT. 

C  **************************************  ************  *******  ii  44  4** 

QLET-10.  **.  01 
DO  t-1.  300 
OPLET-GLET**! 

BPLET*>A<27,  1.  I  ) 

MR ITE (30.  666)  SPLET,  OPLET 
END  DO 

666  FORMAT  (2X.E12.  3.2X.E12.  S) 

END 
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c 

c 
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PROGRAM  UPSET 

******  ******  v*  »«i*******  *******  *******  ****  *4.***************44*** 

***  THIS  PROGRAM  CONUERTS  A  FLUX  VS.  LET  SPECTRUM  INTO  AN 
**#  UPSET /K  IlOI3  IT /DAY  VS.  CRITICAL  CHARGE  PLOT. 

***  REQUIRED  INPUT 

***  X, Y.  Z  »  DIMENSIONS  OF  SENSITIVE  VOLUME 
*«#  REOUIRED  FIl.ES. 

*#*  FOROOO.  DAT  ■  LET  SPECTRUM  FROM  LETINTT 

***  F0R040.DAT  -  OUTPUT  UPSET  PP013  AD  I  LIT  1 ES 

**************  Vi.-  ***************  ******************************* 

DIMENSION  FLUX ( DOG ) . PCI  500  , £R <  500 > . DL ( 500 ) ,  DD ( 500 ) 

REAL  uSOO.LrllN 

osi-a. oat 

***  ACCEPT  DIMENSIONS  OF  SENSITIVE  VOLUME  (MICRONS) 

******  ********  ************************************************ 

ACCEPT  »,.X,  V.  2 

************ ********************************** **************** 

***  COMPUTE  MEAN  PROJECTED  SURFACE  AREA.  CONVERT  TO  PROPER 
***  UNITS. 

********44***************** *********  *************************** 

8-12.  *X*Y+2.  *Y*Z+2.  *2*X1*1.  E-12  !  M**2 

X-X*.  0001 *DS I  !  QRAMS/CM#*2 

V*V*.  0001*DSI  !  QRAM9/CM**2 

2-2*.  0001 *DS I  !  QRAMS/CM**2 

P MAX -SORT ( X*X*Y*Y+2*Z )  !  MAX.  PATHLENGTH  <GM/CM**2) 

**************************************************************** 
***  READ  LET  SPECTRUM 

*■******■**#■****4.-*  *********************************************** 

DO  U-l.  500 

READ  f  30.  400 )  FLUXtJ)-L(J) 

IF  (FLUX(J).  LT.  0.  )  FLUX  ( J  )  ”0. 

FORMAT ( 2 ( 2X.  E12.  S)  ) 

END  DO 

***************•*********#****#******************.************** 

***  COMPUTE  DIFFERENTIALS 

************************************************************** 
01-10.  #*.  01 
02*1 .  -1.  /01 
DO  J-l,  500 

DL  <  J  )  4.02*01  ** J/L ( J )  **2 
END  DO 

******  *44»***»  *****»*****************».**»*  **********4.  *********** 

***  COMPUTE  PATHLENOTH  TERMS 

******  *******  »*«-*****•** ******  »***»***4.*****.»*****4|.** *****»*■*»#* 

DO  J-l.  500 

RAT-10.  **(  Ot«FLOAT( l-J) )  '  RATIO  OF  TWO  STOPPING  POWERS 
DD  (  J  )  — DI FPLD  (  PMAX *RAT,  X,  Y.  2) 

END  DO 

****************************************************4.*****-*** 

»** .  COMPUTE  UPSETS 

**************  4.-********************************************** 

DO  J-l,  500 

QCRIT=PMAX*L'  Jl/22.  S 
PC i J) -OCR  IT 
FACT-22.  5*jCH  IT 
SUM-0 
DO  I  -J,  500 

SUM— SUM+DLI I)*FLUX<I )*DD( i-J+l ) 

END  DO 

ER  (  J  )  — SUM*F  ACT *S*.  25*SS400.  *1024 
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END  DO 

C  *###■*#*•*»*##  »*«■**  #***##**■*  »*******************  «t*****#*****-u-**  *■« 

C  *#*  OUTPUT  UPSET  RATES 

C  *#*■■»*##***##»«■****#«(.**•***♦'■**•»#******■#  #**#****##*#4##*##** **+#4 

DO  J-l ,  SOO 

WRITEUO.  400)  ER  <  J )  i  PC  (  J  > 

END  00 
END 


FUNCTION  D I FPLD ( Si  L.  Ul,  H) 

REAL  L 

AP-3.  *  (H#W+HitL+L»U) 

DIFPLD»(0(9i  L,  W.  H)+0(S.  Wj  L.  H)+a<S,  L.  H>  U)+Q(S.  W.  H.  L)  + 
*  Q<8.  H.  W,L>+Q<3.  H.  L.  W>  >/(3.  14t4*AP) 

END 


FUNCTION  0(9.  X,  V,  Z) 

REAL  KSQ 
KSO»X*X+V#V 
TSQ»X*X-t-Z*Z 
T »3QRT ( TSQ  > 

RSQ-KSG+Z  4Z 
R-SQRT(RSO) 

U-12.  *X*Y#Z*Z 

P9'3«S*S-Z*Z 

QSC5«9*S-X*X-ZJ*Z 

IF  <  <a.  OE.  0.  0).  AND.  <S.  LT.  Z )  >  SOTO  10 
IF((S.  OE.  Z)  .  AND,  (9.  LT.  T)  )  00  TO  SO 
IF<<9.  OE.  T)  .  AND,  (S.  LE.  R>  >  00  TO  30 
0-0.  0 
RETURN 

10  Q-8.  *Y*Y*Z/KSQ-S*<3.  *X*Y/'(R#T!  >4  ■ 

RETURN 

20  0-S<*(3,  *Y/S0RT  !  KSO )  )**2-S*<3.  *X*Y/  ( T*R  >  )  **S! 

*  ~X*<  SQRKPSQ  ) /S  >  *  <9,  +4  *Z»Z/'(S#S)> 

*  +<V*ATAN<Y/'X)-<YfrZ#Z/SURT<KSC5>  )  **2  >  /  ( 3*S*S ) 

RETURN 

30  0— S*(  3,  *X«Z/<R"-SQRT(KSQ)  )  )*-2 

*  +(  X*X*Z*Z*(  Z*Z/KSQ--3.  )+V*ATAN(  Y/X  )  )  X  <S*S*S> 

*  +Y*Z»Z*  (SORT  (CJ5G)  /  9  )*<9.  /TSO+4.  /(S*9>  > 

*  _(V/(3*S*S) ) »AC03(X/SQRT(PSQ) > 

END 
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!  ZCON  C CM  IS  A  COMMAND  FILS  GENERATING  PROGRAM  DESIGNED 

!  FOR  USE  ON  THE  VAX  MINICOMPUTER.  THIS  PROGRAM  WAS  WRITTEN 

!  IN  JANUARY,  1F83  BY  JOHN  R.  LETAW,  SEVERN  COMMUNICATIONS 

1  CQRP  ,  SEVERN*  PARK,  MD 

i 

!  OREETING 

WRITE  SYS4QU  TPIJT  "THIS  COMMAND  FILE  CREATES  A  COMMAND  FILE  WHICH", 
"  WILL  EXECUTE  PROPAGATION  IN" 

WRITE  SYS4GU TPUT  "AIR  PROGRAMS,  AI.l  INPUT  PARAMETERS",  - 
"  ARE  DETERMINED  BY  YOUR  ANSWERS" 

WRITE  SYStOUTPUT  "TO  THE  PROMPTS.  PRESERVE  THIS  SESSION  FOR 
"FUTURE  REFERENCE" 

WRITE  SYStOUTPUT  "  " 

!  COMMAND  FILE  NAME 

INQUIRE  COMPILE  "ENTER  COMMAND  FILE  NAME" 

OPEN/WRITE  OUTF  'COMPILE' 

WRITE  SYS4QUTPWT  11  “ 

!  COMMAND  FILE  DEFAULTS 

WRITE  OUTF  "4  SET  DEFAULT  CLETAW  TEMPI" 

WRITE  OUTF  "4  ASSIGN  1NL:  SYS4PRINT" 

!  INPUT  FLUX 

WRITE  9YG4GUTPUT  "CALCULATE  INPUT  FLUX" 

INQUIRE  YEAR  "ENTER  YEAR  ( 1973.  I44-MIN,  i  9SG,  398-MAX  )  " 

INQUIRE  WEATHER  "ENTER  WEATHER  U-OALACTIC  COSMIC  RAYS)" 

ORBITER-O 

WRITE  OUTF  "4  DELETE  FLUXE3.DAT,*" 

WRITE  OUTF  "4  FOR  FLUX" 

WRITE  OUTF  "4  LIN  FLUX" 

WRITE  OUTF  "4  RUN  FLUX" 

WRITE  OUTF  YEAR,"  ".WEATHER,"  ".ORDITER 

WRITE  OUTF  "4  DELETE  FLUX.  OBJi  *,  FLUX.  EXE)  *,  FLUX.  LIS*  *,  FLUX,  MAP]  #" 

WRITE  SYS40UTPUT  "OUTF  I  LEO.  DAT  CONTAINS  THE  INPUT  FLUX,  " 

PCDUNT-0  1  INITIALIZE  NUMBER  OF  PROPAGATIONS  . 

SUMQRAMS»0  !  INITIALIZE  TOTAL  GfiAMMAGE 

QCOUNT »0  1  INITIALIZE  NUMBER  OF  POST-PROCESSED  FILES 

!  ASSIGN  PROPAGATION  MEDIUM  TO  AIR 

WRITE  OUTF  "4  DELETE  PROP.  *i  *,  RANGIN.  DAT*  #" 

WRITE  OUTF  "4  DELETE  PBWIN.  DAT)  *,  0.  DATi  *" 

WRITE  OUTF  "4  COPY  APRQP.  FQRj  t  PROP  FOR)  1" 

WRITE  OUTF  "4  COPY  ARANGIN,  DAT  I  I  RANQIN.  DAT,  1" 

WRITE  OUTF  "4  COPY  APOWIN.  DAT,  I  POWIN.  DAT,  1 " 

WRITE  OUTF  "4  COPY  AQ.  DAT,  l  Q.  DAT,  1" 

WRITE  OUTF  "4  FOR  PROP" 

WRITE  OUTF  "4  LIN  PROP" 

WRITE  OUTF  "4  FOR  CRAM" 

WRITE  OUTF  "4  FOR  CIJBSPL" 

WRITE  OUTF  "4  LIN  GRAM, CUBSPL " 

WRITE  OUTF  "4  FOR  REFORM" 

WRITE  OUTF  "4  LIN  REFORM" 

PM: 

PCOUNT =P COUNT* I 
!  DETERMINE  GRAHMACE 

WRITE  SYS40UTPUT  "  " 

INQUIRE  GRAMS  "ENTER  NUMBER  OF  GRAMS  Or  AIR  C-D" 
SUMQRAM3-SUMGRAMS+GRAMS 
!  PROP  EXECUTION 

WRITE  OUTF  "4  RUN  PROP" 

WRITE  OUTF  GRAMS 

WRITE  OUTF  "4  DELETE  PROP,  LIS,  *,  PROP.  MAP,  *" 

1  GRAM  EXECUTION 

WRITE  OUTF  "4  DELETE  F0R020.DAT,*" 
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WRITE  OUTF  "3  RUN  ORAM" 

WRITE  OUTF  "*  DELETE  ORAM  LISi  *,  ORAM.  MAP)  *" 

!  ESTABLISH  OUTPUT  FILE 

CUTFILE  -  "OUTF  ILE"-i-FtSTR  INGfPCOUNT )  +  ".  DAT" 

WRITE  OUTF  "1  COPY  FDR020.DAT  ".OUTFILE 
WRITE  SYStOUTPUT  11  11 

WRITE  SYStOUTPUT  OUTFILE."  CONTAINS  THE  FLUX  AFTER  ".SUMQRAMS, 
•'  ORAMS.  " 

!  BRANCH  FOR  FURTHER  PROPAOATIDNS 

INQUIRE  MORE  "DO  YOU  WANT  TO  CONTINUE  PROPAGATION?  <  Y/N ) " 

IF  MORE.  EQS.  "N"  THEN  GOTO  POST 


WRITE  OUTF  "t  DELETE  FL'JXES.  DAT.  «■" 

WRITE  OUTF  "»  RUN  REFORM" 

WRITE  OUTF  n*  DELETE  REFORM.  MAP)  «.  REFORM.  LISi  *" 

GOTO  PM 

!  POST-PROCESSING 

POST: 

WRITE  SYStOUTPUT  "  " 

INQUIRE  MORE  "DO  YOU  WANT  TO  DO  ANY  POST-PROCESSING?  (Y/N)1 
IF  MORE,  EQS.  "N"  THEN  GOTO  END 
I  POST-PROCESSING  OPTION  SELECTION 

WRITE  SYStOUTPUT  "  " 

WRITE  SYStOUTPUT  ~ 

"AT  THIS  POINT  THERE  ARE  SEVERAL  POST-PROCESSING  OPTIONS" 


WRITE  SYStOUTPUT  "  1) 

WRITE  SYStOUTPUT  11  i!> 

WRITE  SYStOUTPUT  "  3) 

WRITE  SYStOUTPUT  "  4) 

WRITE  SYStOUTPUT  "  5) 

INQUIRE  PROCESS  "SELECT  ONE  OPTION" 

!  OPTION  DISTRIBUTION  (1  AND  4) 

IF  PROCESS.  EQS.  1  THEN  GOTO  POST1 
IF  PROCESS.  EQS.  4  THEM  GOTO  P0ST4 
WRITE  OUTF  "t  DELETE  3TP0W.  DAT j *" 

WRITE  OUTF  "t  COPY  SILPOU,  DAT i  i  STPOW.  DAT)  1" 
GOTO  PR0DI3T 

(  OPTION  DISTRIBUTION  <2,  3,  AND  5) 
PRODIST : 

IF  PROCESS.  EQS.  2  THEN  GOTO  POSTS 
IF  PROCESS.  EQS.  3  THEN  GOTO  POSTS 
IF  PROCESS.  EQS.  5  THEN  GOTO  END 
GOTO  POST 
POST1 : 


OMNIDIRECTIONAL  ATMOSPHERIC  INTEGRATION" 
INTEQRAL  LET  SPECTRUM  (TOTAL)" 

INTEGRAL  LET  SPECTRUM  (BY  ELEMENT  GROUP)" 
SOFT  UPSET  RATE" 

EXIT" 


!  ATMOSPHERIC  INTEGRATION  UP  TO  100  QMS  (53.000  FT.  ) 

WRITE  SYStOUTPUT  "  " 

WRITE  SYStOUTPUT  "SEVEN  FILES  CONTAINING  THE  FLUXES  AT  0,  5.  ",  - 

"  15,  30.  50,  75,  AND  100  " 

WRITE  SYStOUTPUT  "GRAHS/CM**2  RESPECTIVELY  ARE  REQUIRED. " 

COUNTER-O 

COUNTS: 

COUNTER »CQUMTER  + 1 

IF  COUNTER.  L ,)  3  THEN  GOTO  XC0UNT2 


INQUIRE  3  "ENTER  FILE1’ 

SS«"AF  ILE"+F»STR  INO  (  COUNTER  )  +>'  DAT" 
WRITE  OUTF  "t  DELETE  ",SS.  "i*" 

WRITE  OUTF  "t  COPY  ".3.  "  " .  SS 


OOTO  COUNTS 


XC0UNT2: 

WRITE  OUTF  " t  FOR  ATMOS" 
WRITE  OUTF  "1  LIN  ATMOS" 
CDEEP-0 
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*  DEPTHS: 

*  INQUIRE  ALTITUDE  "FMTER  ALTITUDE  IN  ATMOSPHERE  IN  KILQFEET" 

*  INQUIRE  CUTOFF  "ENTER  CUTOFF  TIMES  10  IN  GV" 

4  CDESP-CDEEP+1 

*  SS«" ATM'MFPflTR ING  (  ALTITUDE )  +  "  "+F4STR  INS  (  CUTOFF  ) 

*  WRITE  3YS*0UTPUT  S3,"  CONTAINS  THE  FLUX  AT  ".ALTITUDE,"  KILQFEET 

"WITH  A  ",  CUTOFF,  "MO  GV  CUTOFF  " 

*  WRITE  OUTF  "»  RUN  ATMD9" 

*  WRITE  OUTF  ALTITUDE,"  ".CUTOFF 

*  WRITE  OUTF  "*  COPY  FORCxhO  DAT  ",SS 

*  INQUIRE  MORE  "ANOTHER  ALTITUDE  AND  CUTOFF?  ( Y/N  > “ 

*  IF  MORE.  EQS.  "Y"  THEM  GOTO  DEPTHS 

4  WRITE  OUTF  "%  DELETE  AF1LE*.  #i  *,  FDROSO.  DAT,  »" 

*  WRITE  OUTF  "4  DELETE  A TMD3.  EXEi  *.  ATMOS.  ODJj  *,  ATM09.  MAP i  *,  ATMOS.  L IQi  *" 

*  GOTO  POST 
4  POSTS: 

*  !  TOTAL  INTEGRAL  LET  SPECTRUM 

*  WRITE  OUTF  "*  COPY  LGTINTT,  FORi  1  LET.  FOR:  1" 

*  GOTO  CONT 

*  ROOTS’ 

*  !  INTEGRAL  LET  SPECTRUM  QY  ELEMENT  OROUP 

*  WRITE  OUT F  "*  COPY  LETINTG.  FORi  1  LET,  FORM" 

*  GOTO  CONT 

*  P0ST4’ 

*  I  SOFT  UPSET  RATE 

*  WRITE  SYS40UTPUT  11  " 

*  WRITE  SYS40UTPUT  "CONUERT  AN  INTEGRAL  LET  FILE  INTO  SOFT  UPSET  RATE" 

*  WRITE  SYS40UTPUT  "  " 

*  WRITE  OUTF  "4  FOR  UPSET" 

4  WRITE  OUTF  "*  LIN  UPSET" 

*  COUNTERED 

*  COUNT: 

*  COUNTER-COUNTER#! 

4  SS*"EF 1LE  "’"FtSTR  I  NO  <  COUNTER  >  +  ",  DAT" 

»  INQUIRE  S  "INTEGRAL  LET  FILE" 
t  WRITE  SYS*OUTPUT  "  " 

*  WRITE  OUTF  "*  DELETE  F0R030.  DAT  I  *" 

*  WRITE  OUTF  "t  COPY  », 8, "  FQR030.DAT, 1" 

*  WRITE  OUTF  "*  DELETE  F0RO4Q.  DAT:  *“ 

*  WRITE  9YS*0UTPUT  "INPUT  DIMENSIONS  OF  SENSITIVE  REGION  (MICRONS)" 

*  INQUIRE  DIM!  "LENGTH" 

*  INQUIRE  DIMS  "WIDTH  " 

*  INQUIRE  DIM3  "DEPTH  " 

*  WRITE  3YS40UTPUT  "  " 

*  WRITE  OUTF  "S  RUN  UPSET" 

*  WRITE  SYS40UTPUT  "UPSET  RATE  FOR  ",S. "  IS  IN  " ,  S3 

*  WRITE  OUTF  DIM1,"  ".DIMS,"  ",DIM3 

*  WRITE  SYS40UTPUT  "  " 

*  WRITE  OUTF  "4  COPY  FOR040  DAT  ",  S3 

*  WRITE  OUTF  "I  DELETE  FQR040  DAT;  *,  FOR030.  DAT:  *" 

*  WRITE  OUTF  " »  DELETE  UPSET  LI  Si#, UPSET  MAP, 

*  INQUIRE  MORS  "MORE  UPSET  CALCULATIONS?  <Y/N>" 

*  WRITE  3YS40UTPUT  "  " 

4  IF  MORE,  EOS.  "Y"  THEN  GOTO  COUNT 

*  WRITE  OUTF  "»  DELE Tu  UPSET.  EXEi  *,  UPSET.  ODJi  #“ 
t  GOTO  POST 

1  CONT; 

*  WRITE  OUTF  "*  FOR  LET" 

*  WRITE  OUTF  "1  FOR  CUnSPL" 

*  WRITE  OUTF  "4  LIN  LEV, CUGSPL" 

*  PFILES: 
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*  WRITE  SYS*OUTPUT  "  " 

»  INQUIRE  FILENAME  "FILE  TO  DE  PROCESSED" 

*  QCOUNT-QCOUNT+1 

*  WRITE  OUTF  “*  DELETE  FOR020.  DATi *" 

*  WRITE  OUTF  "*  COPY  "■FILENAME,"  F0R020. DAT i 1 " 

*  WRITE  OUTF  "*  RUN  LET" 

*  POSFILE  ■  llPFILE"+F*STR  ING(  QCOUNT)  +" .  DAT" 

*  WRITE  OUTF  "*  COPY  FDR030. DATi  I  ",  POSFILE 

*  WRITE  OUTF  "*  DELETE  FQR030.  DATi #" 

*  WRITE  SYSiOUTPUT  FILENAME,"  WAS  PROCESSED  INTO  ".POSFILE 

*  INQUIRE  MORE  "ANY  MORE  FILES  TO  BE  PROCESSED  THIS  WAY?  <V/N)" 

*  IF  MORE. EQS. "Y"  THEN  OOTO  PFILES 

1  WRITE  OUTF  "*  DELETE  LET.  *i  *,  CUBSPL.  ODJi  #,  CUESPL.  LISi  *" 

»  QOTO  POST 

*  END: 

*  WRITE  OUTF  "*  DELETE  FOR*.  *i  #,  STPOW.  DATi  *" 

*  CLOSE  OUTF 

*  EXIT 
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